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SUMMARY 


Presented  in  this  investigation  is  a  theoretical  analysis  of 
the  thrust  augmentation  characteristics  of  jet  ejectors. 

The  analysis  includes  the  effects  of  flow  compressibility, 
major  flow  losses,  and  forward  speed.  Numerical  rr  ults  are 
presented  in  the  form  of  nomographs  for  a  wide  range  of 
practical  operating  conditions.  These  computations  were 
performed  with  the  aid  of  an  IBM  360  digital  computer.  The 
chaits  can  be  used  to  predict  the  jet  ejector  performance 
and  as  such  represent  an  effective  analytical  tool  for 
preliminary  design  purposes.  The  numerical  results  are  used 
to  determine  the  effects  of  the  more  important  aerodynamic, 
thermodynamic,  and  geometric  parameters  on  jet  ejector 
thrust  augmentation.  A  correlation  of  these  results  with  the 
available  experimental  data  is  also  made. 
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I.  INTRODUCTION 


The  concept  of  generating  thrust  augmentation  by  means  of  a 
jet  ejector  became  an  active  research  subject  several 
decades  ago.  Since  then,  numerous  studies  have  been 
conducted  providing  a  considerable  amount  of  jet  ejector 
data.  Much  of  this  work  was  done,  however,  for  different 
and  diverse  applications,  and  because  of  the  complexity  of 
the  jet  ejector  flow  problems,  the  available  information 
cannot  be  readily  used  for  practical  applications  to  aircraft 
design  and  performance. 

The  main  objective  of  this  program  is  to  evaluate  the  avail¬ 
able  data  as  to  their  practical  applicability  and  to  extend 
and  modify  existing  theories  to  provide  an  analysis  for 
realistic  appraisal  of  the  potential  of  jet  ejectors  for 
achieving  augmented  thrust. 

The  review  of  the  existing  literature  covered  a  major 
portion  of  technical  reports  (total  585)  which  are  listed 
in  Appendix  II.  A  discussion  of  the  more  important  of  these 
investigations  is  presented  in  Section  III. 

The  theoretical  analyses  which  are  formulated  in  Section  IV 
are  performed  for  an  axisymmetric  ejector  and  include  the 
effects  of  compressibility,  flow  losses,  and  ejector  geom¬ 
etry.  These  analyses  utilize  various  simplifying  assumptions, 
among  which  are  the  conditions  that  the  velocity  profile  at 
the  secondary  entrance  is  uniform  and  that  the  mixing  of  the 
primary  and  secondary  flows  is  completed  at  the  exit  of  the 
mixing  chamber.  The  former  assumption  limits  the  appli¬ 
cability  of  this  one-dimensional  theoretical  approach  to 
not  too  large  secondary-to-primary  area  ratios.  A  three- 
dimensional  analysis  is  indeed  a  formidable  task  and  has  not 
been  undertaken  by  any  of  the  investigators  in  the  past.  To 
keep  the  complexity  of  the  analysis  within  the  scope  of  the 
present  program, an  empirical  correction  factor  for  thrust 
augmentation  ratio  is  herein  formulated  to  account  for  the 
effect  of  the  nonuniform  velocity  distribution  at  the 
secondary  entrance.  The  latter  assumption  implies  that  the 
mixing  chamber  length  of  a  given  ejector  configuration  must 
be  sufficiently  long  to  ensure  complete  mixing  of  the  exit 
flows.  No  precise  information  is  available  on  this  subject, 
and  therefore  a  semiempirical  approach  is  herein  utilized 
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for  qualitative  evaluation  of  the  mixing  chamber  length 
required  for  complete  mixing. 

Furthermore,  a  comprehensive  parametric  study  is  conducted 
to  determine  the  effects  of  the  more  important  aerodynamic, 
thermodynamic,  and  geometric  parameters  on  jet  ejector 
performance.  A  discussion  of  these  effects  is  presented  in 
Section  V. 

Section  VI  contains  a  correlation  of  the  theoretical  results 
with  the  available  experimental  data. 

Section  VII  contains  a  summary  of  the  theoretical  results 
for  rapid  predictions  of  jet  ejector  performance.  These 
results  were  obtained  with  the  aid  of  an  IBM  360  digital 
computer  and  are  presented  as  nomographs. 
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II.  JET  EJECTOR  PRINCIPLE 


A  jet  ejector  is  a  device  in  which  a  secondary,  or  driven, 
fluid  is  entrained  by  a  primary  or  actuating  fluid  with 
subsequent  transfer  of  energy  through  turbulent  mixing  in  a 
mixing  chamber.  The  primary  fluid,  which  is  originally  at  a 
higher  stagnation  pressure,  is  discharged  with  a  high 
velocity  into  the  mixing  chamber  of  specific  shape.  Due  to 
viscous  shear,  the  fluid  surrounding  the  primary  flow  is 
brought  into  motion  at  the  entrance  of  the  mixing  chamber. 
This  motion  causes  a  drop  of  static  pressure,  as  a  result  of 
which  the  secondary  fluid,  in  many  cases  ambient  air,  is 
entrained  into  the  mixing  chamber.  The  secondary  flow  thus 
formed  mixes  turbulently  with  the  primary  jet  in  "the  mixing 
chamber  and  energy  transfer  occurs.  The  mixed  flow  then 
proceeds  toward  the  exit  end  of  the  mixing  chamber  and 
finally  discharges  to  some  back  pressure  which  may  be 
atmospheric.  If  a  diffuser  is  attached,  the  mixed  flow 
builds  up  some  static  pressure  before  reaching  the  exit.  As 
a  result  of  the  pumping  action  as  described  above,  the  total 
momentum  of  the  mixed  flow  at  the  ejector  exit  is  increased 
due  to  the  entrainment  of  the  secondary  fluid,  as  compared 
with  the  momentum  of  the  primary  jet  discharged  directly  into 
the  atmosphere.  Jet  thrust  augmentation  is  thus  achieved. 

In  essence,  the  jet  ejector  can  be  considered  as  a  device 
which  converts  a  low  mass  flow  propelled  at  high  velocity  to 
a  high  mass  flow  propelled  at  low  velocity.  If  this  is 
accomplished  at  little  energy  loss,  thrust  augmentation  can 
be  achieved. 

The  efficiency  of  this  conversion  process  depends  on  the 
ejector  configuration,  the  detail  geometry  of  these  configu¬ 
rations,  and  the  thermodynamic  and  aerodynamic  operating 
conditions . 

The  more  important  jet  ejector  configurations  are  as  follows: 

(a)  Single  nozzle  ejector 

(b)  Multiple  nozzle  ejector 

(c)  Annular  nozzle  ejector 


The  prime  parameter  affecting  the  performance  of  these 
configurations  is  the  minimum  mixing  chamber  length  required 
for  complete  mixing  of  the  primary  and  the  secondary  flows. 
No  precise  analytical  methods  are  available  for  determining 
this  parameter;  therefore,  a  semiempirical  method  is  herein 
formulated  for  this  purpose.  An  analysis  of  various  ejector 
types  is  presented  in  Section  IV. 

The  major  geometric  parameters  affecting  thrust  augmentation 
are : 

(a)  Secondary-to-primary  area  ratio  at  the  entrance 
to  the  mixing  chamber. 

(b)  Shape  and  length  of  the  mixing  chamber. 

(c)  Diffuser  exit-to-entrance  area  ratio  and  diffuser 
angle . 

(d)  Primary  nozzle  configuration  and  location. 

(e)  Secondary  entrance  inlet  contour. 

Finally,  the  thermodynamic  and  aerodynamic  parameters  of 
importance  are: 

(a)  Stagnation  properties  of  the  primary  fluid. 

(b)  Stagnation  properties  of  the  secondary  fluid. 

(c)  Atmospheric  conditions  at  the  ejector  exit. 

(d)  Forward  speed  of  the  ejector  system. 
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III. 


REVIEW  OF  THE  AVAILABLE  PERTINENT  LITERATURE 


Presented  in  this  section  is  a  brief  review  of  the  state  of 
the  art  of  jet  ejectors.  The  papers  selected  for  this 
discussion  are  those  which  are  more  directly  related  to  the 
work  under  the  present  program,  and  they  are  herein  arranged 
in  a  chronological  order.  A  more  complete  bibliography  on 
jet  ejectors  is  presented  in  Appendix  II. 

The  analysis  performed  by  Roy  (Reference  1)  deals  with 
compressible  fluid  flow  under  practical  operating  conditions 
of  jet  ejeccors.  Entrance  and  diffuser  losses  are  accounted 
for  by  assuming  the  flow  processes  as  polytropic.  The 
friction  loss  at  the  mixing  chamber  wall  is  also  considered. 
This  analysis  results  in  a  system  of  24  flow  equations  for 
which  no  close-form  or  digital  computer  solutions  are 
presented. 

Morrisson  (Reference  2)  performs  an  incompressible  analysis 
in  which  it  is  assumed  that  the  mass  flow  rate  of  the 
primary  jet  exhausted  to  the  ambient  is  the  same  as  that 
discharged  into  the  mixing  chamber.  This  assumption  may  not 
be  valid  since  the  pressure  reduction  created  in  the  mixing 
chamber  results  in  an  increase  of  mass  flow  rate  of  the 
primary  jet  as  compared  to  that  when  the  primary  jet  is 
discharged  into  the  atmosphere.  The  analysis  presents 
charts  of  thrust  augmentation  ratio  versus  the  mass  flow 
ratio.  However,  these  charts  cannot  readily  be  applied  for 
determining  jet  ejector  performance,  since  the  mass  flow 
rate  is  unknown. 

Sargent  (References  3  and  4)  performs  idealized  analyses  of 
the  performance  of  jet  ejectors  for  both  static  and  forward 
speed  conditions.  In  these  analyses,  the  flows  are  assumed 
to  be  incompressible  and  the  solution  for  the  flow  equations 
is  made  by  a  trial  and  error  method  with  the  area  ratio  and 
mass  flow  ratio  as  variables.  These  analyses  rely  on  the 
assumption  that  the  mass  flow  rate  of  the  primary  jet 
exhausted  to  the  ambient  is  equal  to  chat  of  the  primary  jet 
inside  the  ejector. 

The  analysis  of  McClintock  et  al  (Reference  5)  results  in  an 
e  Ration  for  thrust  augmentation  as  a  function  of  the  mass 
flow  ratio  which  is  not  determined.  In  spite  of  the  fact 
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that  the  analysis  is  performed  for  incompressible  fluids,  an 
attempt  is  made  to  account  for  different  densities  of  primary 
and  secondary  flows.  This  appears  to  be  contradictory,  since 
an  incompressible  analysis  is  only  valid  for  primary  and 
secondary  fluids  of  nearly  equal  densities. 

Ellerbrock's  analysis  (Reference  6),  which  includes  the  loss 
due  to  wall  friction,  treats  compressible  flow  for  constant 
area  mixing.  However,  in  obtaining  solutions  to  the  flow 
equations,  assumed  values  of  the  pressure  ratio  at  the 
ejector  entrance  plane  to  that  of  the  ambient  are  introduced. 

The  analysis  performed  by  von  Karman  (Reference  7)  presents 
the  basic  information  on  cylindrical  jet  ejectors  operating 
under  idealized,  incompressible  flow  conditions.  The 
expressioi.  for  the  nondimensionalized  velocity  of  the  mixed 
flow  at  the  exit  appears  to  be  erroneous,  although  the 
relationship  xor  thrust  augmentation  ratio  finally  obtained 
is  correct.  in  this  analysis,  an  attempt  is  made  to  account 
for  the  effect  of  the  nonuniform  velocity  profile  at  the 
secondary  entrance  on  jet  ejector  thrust  augmentation.  How¬ 
ever,  the  analysis  utilizes  the  assumption  that  the  mass  flow 
rate  for  the  nonuniform  velocity  profile  is  the  same  as  that 
in  the  case  of  uniform  velocity  distribution..  As  the  result 
of  this  assumption,  an  increase  rather  than  a  reduction  in 
thrust  augmentation  is  achieved. 

Szczeniowski. ,  in  his  incompressible  analysis  (Reference  d), 
considers  a  hypothetical  approach  in  which  momentum  and 
energv  are  transferred  between  the  primary  and  secondary 
streams  while  the  two  flows  remain  separate  at  the  exit  of 
the  ejector.  The  justification  of  this  approach  is  rather 
cl  i  f  f  icult  . 

Sandv '-s  et.  al  (Reference  9)  compute  the  ejector  thrust  by 
integ.ation  of  surface  pressures;  however,  the  momentum  of 
the  flow  from  the  exit  of  the  ejector  is  completely  ignored. 
This  approach  does  not  yield  the  total  thrust  augmentation. 

The  analysis  performed  by  Bert  in  et  al  (Reference  10) 
utilizes  the  usual  one-dimensional  flow  approach  but  yields 
no  close- form  solution  for  the  jet  ejector  performance.  The 
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performance  charts  presented  in  this  reference  have  limited 
practical  application  because  ehe  mass  flow  ratio  whicn  is 
used  as  an  independent  variable  cannot  be  determined. 

Chisholm's  work  (Reference  11)  covers  compressible  and 
incompressible  flows  with  constant  area  and  constant  pressure 
mixing.  Although  the  ejector  flow  equations  are  formulated, 
no  solution  has  been  attempted  for  determining  tnc  thrust 
augmentation  ratio. 

Reid's  investigation  (Reference  12)  is  chiefly  aimed  at  the 
reduction  of  jet  noise  by  means  of  an  ejecto*-;  however,  an 
expression  of  a  so-called  thrust  parameter  is  presented.  A 
simplified  theoretical  analysis  foi  constant  i  mixing  with 
assumed  common  stagnation  temperature  for  the  primary  and 
secondary  streams  is  presented.  The  analysis  indicates  that 
although  the  system  of  flow  equations  formulated  can  be 
solved  in  principle,  the  numerical  solutions  are  very 
difficult  to  obtain. 

The  analysis  presented  in  Reference  13  is  for  a  constant 
pressure  mixing  jet  ejector  with  diffuser.  The  primary  and 
secondary  flows  are  assumed  to  have  the  same  stagnation 
temperature.  As  pointed  out  in  this  reference,  the  system 
of  equations  can  be  solved  when  the  mixing  pressure  is  known. 
However,  since  the  mixing  pressure  is  one  of  the  parameters 
to  be  determined,  this  analysis  has  no  direct  practical 
application. 

Storkebaum  in  his  work  (Reference  iw)  deals  with  both 
incompressible  and  compressible  flow  analyses.  For  tne 
incompressible  case,  the  idealized  thrust  augment  at  -.on  ratio 
is  presented  as  a  function  of  the  secondary- to-primary  area 
ratio.  The  effect  of  the  losses  at  the  entrance  to  the 
mixing  chamber  is  uiso  uuied.  In  the  compressible  analysis, 
the  equations  for  thrust  augmentation  are  formulated; 
however,  the  procedure  of  solving  the  equations  is  not 
presented . 

Wan  performs  an  incompressible  analysis:  (Reference  15) 
similar  to  that  presented  by  von  Karma n  (reference  7),  with 
the  exception  chat  the  mixed  flew  velocity  at  the  ejector 
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exit  is  considered  to  be  nonuniform.  However,  the  validity 
of  this  analysis  appears  to  be  doubtful  because,  despite  the 
nonuniformity  of  the  exit  velocity,  the  static  pressure  at 
the  exit  is  assumed  to  be  uniform. 

Payne  (Reference  16)  conducts  a  theoretical  analysis  of  a 
jet  ejector  with  constant  pressure  mixing  treating  flows  as 
incompressible.  This  analysis  indicates  that  the  optimum 
thrust  augmentation  is  primarily  dependent  on  the  diffuser 
efficiency  and  that  augmentation  ratios  as  high  as  4.0  or 
more  are  possible  with  high  diffuser  efficiencies.  This 
conclusion  is  based  on  an  infinite  secondary-to-primary  mass 
flow  ratio  which  cannot  be  achieved  in  practice.  Also,  in 
performing  the  differentiation  process  to  obtain  optimum 
thrust  augmentation,  the  analysis  uses  the  pressure  in  the 
mixing  chamber  as  the  independent  variable  and  considers  the 
mass  flow  ratio  as  a  constant.  This  appears  to  be  incorrect, 
since  the  mass  flow  ratio  appearing  in  the  basic  flow 
equation  is  a  function  of  the  pressure  in  the  mixing  chamber 
and  cannot  be  treated  as  a  constant.  Furthermore,  in  the 
breakdown  of  the  Location  of  thrust  increase,  the  analysis 
indicates  that  the  largest  contribution  to  the  total  thrust 
augmentation  is  due  to  the  bellmouth  lip  intake  rather  than 
the  change  of  momentum  of  incoming  and  outgoing  flows.  No 
experimental  data  are  presented  to  justify  this  conclusion. 

Sandover  (Reference  17),  utilizing  the  work  of  Payne 
(References  16  and  18),  performs  an  analysis  also  based  on 
one-dimensional  incompressible  fluid  flow.  This  analysis 
commences  with  an  assumption  that  the  constant  pressure  in 
the  mixing  chamber  can  be  achieved  with  a  uniform  cross 
section  of  the  chamber;  then  an  attempt  is  made  to  determine 
the  optimum  exit-to-entrance  area  ratio.  Furthermore,  the 
analysis  indicates  that  an  increase  in  thrust  augmentation 
can  be  achieved  by  increasing  the  temperature  of  the  primary 
jet  flow.  This  is  considered  to  be  incorrect  for  both  the 
incompressible  or  compressible  flow  analyses.  In  the  former 
case,  temperature  does  not  affect  tae  thrust;  in  the  latter 
case,  a^  indicated  in  the  present  report,  the  increase  of 
the  primary  jet  temperature  results  in  a  decrease  ratner  than 
an  increase  of  the  thrust  augmentation  ratio.  This  result 
has  been  verified  by  available  experimental  data. 
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From  Che  above  review,  it  is  seen  that  little  usable  informa¬ 
tion  on  jet  ejector  performance  is  available  in  the  existing 
technical  literature.  Despite  the  fact  that  in  numerous 
investigations,  equations  representing  the  ejector  flows  are 
formulated,  no  explicit  solutions  are  presented  for  the 
practical  evaluation  of  jet  ejector  performance.  In  the 
cases  where  some  numerical  solutions  are  presented,  they  are 
generally  expressed  as  functions  of  such  parameters  as  the 
mass  flow  ratio,  the  pressure  at  the  mixing  chamber  entrance, 
primary  jet  Mach  number,  etc.,  none  of  which  can  be  readily 
determined. 


/■* 
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IV.  THEORETICAL  ANALYSES 


As  pointed  out  in  the  previous  section,  the  presently 
available  literature  does  not  provide  adequate  information 
on  the  thrust  augmentation  characteristics  of  jet  ejectors. 
Existing  analyses  were  therefore  modified  and  extended  and 
practical  methods  were  formulated  by  means  of  which  the 
potential  of  jet  ejectors  can  be  effectively  evaluated. 

Because  of  the  complexity  of  the  problem,  the  following 
simplifying  assumptions  are  made: 

(a)  The  ejector  geometry  is  axi symmetrical . 

(b)  The  vt’ocity  and  static  pressure  of  the  primary 
and  secondary  flow  are  uniform  at  the  entrance 
plane  to  the  mixing  chamber. 

(c)  The  static  pressures  of  the  primary  and  secondary 
flows  at  the  entrance  to  the  mixing  chamber  are 
equal . 

(d)  The  velocity  and  static  pressure  of  the  mixed  flow 
are  uniform  at  the  exit  of  the  mixing  chamber,  and 
in  the  case  of  an  ejector  with  diffuser  are 
uniform  for  any  cross  section  of  the  diffuser. 

(e)  Both  the  primary  and  secondary  fluids  are  perfect 
gases  having  the  same  specific  heat. 

(f)  The  pressure  at  the  exit  of  the  mixing  chamber,  or 
in  the  case  of  a  diffuser,  at  the  exit  of  the 
diffuser,  is  equal  to  the  ambient  pressure. 

(g)  The  stagnation  conditions  of  the  primary  flow  are 
unaffected  by  the  presence  of  the  ejector,  and  the 
losses  in  the  primary  flow  up  to  the  exit  of  the 
nozzle  are  neglected. 

(h)  No  heat  losses  occur  at  the  mixing  chamber  and 
diffuser  walls. 

(i)  The  motion  of  the  ejector  system,  if  any,  is  in 
the  same  direction  as  that  of  r'.e  primary  jet  flow. 
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With  these  simplifying  assumptions,  the  ejector  flow  field  is 
reduced  to  a  one-dimensional  fluid  flow  pioblem  which  can  be 
treated  analytically  utilizing  the  following  principles  of 
fluid  dynamics: 

(a)  Conservation  of  mass. 

(b)  Conservation  of  energy. 

(c)  Newton’s  second  law. 

(d)  Equation  of  state  for  perfect  gases. 

(e)  Thermodynamic  processes  of  the  flows. 

The  formulation  of  the  jet  ejector  flow  equations  is  herein 
accomplished  utilizing  Figure  1.  As  indicated  in  this 
figure,  various  stations  of  jet  ejectors  are  designated  by 
numbers  (0),  (1),  (2),  and  (3)  which  are  assigned  to 
represent  the  respective  stations  of  the  jet  ejector  flows 
as  follows: 


Station 

(0) 

Upstream  of  the  entrance  of  the  mixing 
chamber. 

Station 

(1) 

Plane  at  the  entrance  of  the  mixing 
chamber,  which  is  also  the  exit  plane 
of  the  primary  nozzle. 

Station 

(2) 

Plane  at  the  exit  of  the  mixing  chamber, 
which  is  also  the  plane  at  the  entrance 
of  the  diffuser. 

Station 

(3) 

Plane  at  the  exit  of  the  diffuse. 

Utilizing  the  nomenclature  of  Figure  1  and  the  assumptions 
and  conditions  described  above,  the  following  jet  ejector 
flow  equations  are  obtained: 
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Primary  Flow 


Secondary  Flow 
Mixed  Flow 


Figure  1.  Schematic  Representation 
of  Jet  Ejector  Configuration. 
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Conservation  of  Energy  for  the  Primary  Flow  up  to  the 
Nozzle  Exit 


CP  T0p  J 


_r__P <_  %2 

r*1  />ip  +  2g 


(6) 


Conservation  o£  Energy  for  the  Secondary  Flow  up  to 
the  Entrance  to  the  Mixing  Chamber 


cp  Tq  J 


r_  < 

X'1  Pis  2g 


(7) 


Isentropic  Process  of  the  Primary  Flow 

- (8! 

p0p  Vr1  ^ipCpTopJ/ 

Irreversible  Adiabatic  Process  for  the  Secondary  Flow 
up  to  the  Entrance  to  the  Mixing  Chamber 


I  y  P, 

—  _ - ! - +  I  ;  - 

Ve  y1  PiscpTqJ  \  Ve 


2g  CpigJ 


y-1  (9) 


Irreversible  Adiabatic  Process  in  the  Dit f u s e r 


_P_2 

Pn 


_ _ Pz_  ,  V2 

P2  l  y-\  pz  2  g 


,  r  .  P *2_  ,  V2' 

p*\y-\  p2 


2g 
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The  unknown  quantities  in  the  above  equations  are:  P(  ,  P?  , 


equations  with  10  unknowns.  As  can  be  noted,  no  close-form 
solution  for  these  equations  is  possible,  and  even  an 
iterative  solution  with  the  aid  of  a  digital  computer  is 
indeed  a  formidable  task. 

Hence,  in  order  to  obtain  a  practical  evaluation  of  jet 
ejector  performance,  further  simplifications  are  necessary 
in  the  analysis. 

The  present  approach  consists  of  first  formulating  a  very 
simplified  analysis,  defined  herein  as  the  idealized 
analysis,  which  enables  rapid  prediction  of  maximum  thrust 
augmentation  and  the  entrainment  ratio  of  jet  ejectors,  for 
the  most  idealized  flow  conditions.  The  idealized  analysis 
is  herein  formulated  considering  perfect  fluid  flow  with  no 
losses  and  no  forward  speed  effects.  This  analysis  is 
performed  for  two  types  of  jet  ejectors  with  constant  area 
mixing  and  constant  pressure  mixing. 

As  a  second  step,  a  practical  analysis  is  formulated.  This 
analysis  consists  of  two  parts: 

(a)  The  effects  of  major  flow  losses,  forward  speed, 
and  diffuser  are  considered  assuming  incompressible 
fluid  flow. 

(b)  The  effect  of  flow  compressibility  is  determined 
on  the  basis  of  no  flow  losses,  no  diffuser,  and 
no  forward  speed. 

This  approach  makes  it  possible  to  perform  with  reasonable 
computational  effort  a  comprehensive  parametric  study  of  the 
relative  importance  of  the  individual  flow  parameters  on  jet 
ejector  performance  and  also  provides  a  practical  method  for 
the  evaluation  of  various  jet:  ejector  systems. 

A.  IDEALIZED  ANALYSIS 


The  idealized  analysis,  which  in  part  is  also  presented  in 
Reference  7,  utilizes  the  following  additional  simplifying 
assumptions : 
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(a)  The  flows  are  incompressible,  which  implies  that 
the  fluid  density  is  uniform  and  constant 
throughout  for  all  the  flows  and  is  unaffected  by 
temperature  and  pressure. 

(b)  All  flow  losses  are  neglected. 

(c)  The  ejector  system  has  no  forward  speed. 

These  assumptions  apply  to  both  constant  area  mixing  and 
constant  pressure  mixing. 

1 .  Constant  Area  Mixing 

The  thrust  augmentation  ratio  is  defined  as  the  ratio  of 
the  thrust  produced  with  an  ejector  on  F3  to  that  produced 
by  the  primary  jet  when  it  is  discharged  directly  into  the 
atmosphere  F Qp. 

Thus , 


(11) 


where  the  thrust  produced  by  the  ejector  is  given  by 


F 


3 


P_ 

q 


A  3 


(12) 


and  the  thrust  produced  by  the  primary  jet  discharged 
directly  into  the  atmosphere  is 


t  A 

q  Aip 


Vor 


(13) 
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Substituting  equations  (12)  and  (13)  into  equation  (11), 
there  follows 


<t>  - 


a  3  Vj 

Aip  Vap 


(1*) 


Applying  Bernoulli's  equation  for  the  primary  flow 
discharged  into  the  atmosphere  and  the  primary  flow 
discharged  into  the  ejector,  and  assuming  that  the 
stagnation  pressure  P0p  is  not  affected  by  the  presence  of 
the  ejector,  there  results 


POp:PQ 


van-'  pm  + 


(15) 


Also  applying  Bernoulli's  equation  between  stations  (0)  and 
(1)  for  the  secondary  flow, 


(16) 


From  equations  (15)  and  (16)  there  results 


P|  +  c  V,  +  ~  p  V0  -  P,  -f  ——  p  V, 
s  cq'  s  2q  r  aP  p  2  g  ■  1  p 


(17) 


Using  the  assumption  that  Plp  *  P,$  ,  equation  (17)  reduces  t 


2  2  2 
Vap  ;  %  ~vis 
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or 


Substituting  equation 


(18)  into  (14)  yields 


(18) 


4>- 


(19) 


The  analysis  will  now  proceed  to  develop  the  relationships 
for  V i $ / V , p  and  V3/V|_  in  terms  of  the  area  ratios 
aE:  A,s/Alp  'uul  a  o'  A3/A2  =  A3^A(p  +Ais) 

dsing  the  continuity  equation  between  stations  (1)  and  (3), 
there  results 


A  V.  +A,  V.  =A,V 

■p  'p  's  's 


3  ’  3 


Kquation  130)  reduces  to 


/  V  3  \ 

toc  +  iloo^-j-l 


at 


(30) 


( 3 1  ) 


Applying  Bernoulli’s  equation  between  stations  (2)  and  (3), 
vising  P3  PQ  ,  there  follows 


(23) 
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Substituting  equation  (16)  into  equation  (22)  yields 


p^  +  ^C+-^’v>2--27^2  <23 

Also,  from  momentum  considerations  between  stations  (1)  and 
(2),  there  results 


p  +  ^'s  ^'S  ^2^2  :  g  ^2^: 


—  A,  V,  2 — —  A|  V,2  (24) 
q  'p  'p  q  's  v,s 


Substituting  equation  (23)  into  (24)  and  using  P,n=P 
yields 


p  '  '$ 


i  V  WA*,fV V  2 kpV’~ ^pV‘} - Ta^- |  VV  f‘ W* 


{ 2  0 , 


Equation  (25)  can  be  simplified  as  follows 


\cic  +0 


-2 


(26) 


From  continuity  equation  between  stations  (2)  and  (3), 


V2  =aDV3 


(27) 


Squaring  equation  (27)  and  substituting  into  equation  (26) 
vie  1 ds 

2 

2  (aE-H)(a0+D 


at-l 
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Squaring  equation  (21)  and  equating  it  to  equation  (28), 
there  results  the  following  quadratic  equation  for  Vj/V^  : 


/\M_  2(aEH>aD 

Up/  a^-ao2- 


Vlp  /  aE+aD 


(29) 


Solving  equation  (29),  the  relationship  for  V3/V,  is  as 
follows: 


-(aE-l)aD±  aE  Ua0  +  2aE- 


a 


2 

D 


(30) 


It  can  be  noted  that  the  velocity  ratio  V3/V1  p  must  always 
be  greater  than  zero  for  all  positive  values  of  the  area 
ratios  aE  and  aD.  Since  a0  >  1,  there  fellows  that  for 
aE  >  0,  the  positive  sign  in  front  of  the  second  terms 

must  be  used.  Hence, 


/  2 

-(aE-l)ao  +  aE  V  oto  +  2aE- 

2,  2 

ceE  +  aD 


(31) 


Substituting  equation  (31)  into  equation  (21)  yields 

r  _ "i 


(a 


;+l)aD  -(aE-!)a0  +  aEVa^  '+  2 aE-l  J  -(aE2+  ao) 


aE(aE2+aQ) 


(32) 


Finally,  substituting  equations  (31)  and  (32)  into  elation 
(19),  the  thrust  augmentation  ratio  <p  can  be  analytically 
expressed  as  follows: 


4>z 


aE  (aE+l)a0 


(aE~l)  a0 


'  2  ~ 

+  aEV  a0  +  2aE 


t 

-I  J 


2  2  2  2  r  r  /~Z  -1  2  2^2 

aE  (aE  +aD)-  |(aE  +  l)a0[-(c:E-l)  aD+aEv  aD-r2aE-l  J  -(aE  +  ao)j 


(33) 
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Similarly,  the  mass  entrainment  ratio  w  is  given  by 


W 


\% 


[/  2  ^  2  2 
~  1 )  +  &£  V  &q"^2Q£~Ij  “  ( Gt  £  CE  q  ) 


aE2+a02 


(34) 


For  a  special  case  of  constant  area  mixing  condition  with  no 
diffuser,  i.e.,aD  -=  1.0,  the  thrust  augmentation  ratio  <f> 
reduces  to 


(QE  + 


laE+  l)[-laE-l)  +  a 

1  2  q  +  (a 


»y^r 


and  the  mass  entrainment  ratio  becomes 


(35) 


W  = 


(a. 


->[■ 


(aE-l)  -f  aE  y  2aEj  ~laE  +  l) 


aE+  I 


(36) 


It  can  be  seen  from  equation  (35)  that  <p  approaches  a  limit 
of  2.0  as  aE  tends  to  infinity.  This  obviously  is  not  true 
even  under  idealized  conditions,  since  an  ejector  with 
infinite  secondary-to-primary  area  ratio  is  equivalent  to  a 
free  jet  discharged  into  the  atmosphere,  in  which  case  no 
thrust  augmentation  exists  and  hence  <fi  =  1.0.  This 
discrepancy  is  attributed  to  the  unrealistic  assumption  that 
the  secondary  velocity  at  the  entrance  to  the  mixing  chamber 
is  uniform  regardless  of  the  area  ratio.  A  discussion  on 
this  subject,  together  with  the  derivation  of  an  empirical 
correction  factor  for  the  no:  uniform  secondary  entrance 
velocity  profile,  is  presented  later  in  this  report. 

2 .  Constant  Pressure  Mixing 

The  mixing  chamber  of  an  ejector  is  sometimes  of  variable 
cross-sectional  area.  Due  to  the  lack  of  information  on  the 
longitudinal  pressure  variation  inside  the  mixing  chamber, 
it  is,  in  general,  difficult  to  express  the  pressure  force 
on  the  walls  in  the  momentum  equation.  However,  one 
configuration  exists  for  which  an  analysis  can  easily  be 
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performed.  This  is  the  so-called  constant  pressure  mixing 
configuration  which  implies  that  the  pressure  inside  the 
mixing  chamber  is  considered  constant  throughout. 

Thus,  for  constant  pressure  mixing,  the  additional  condition 
is  that  the  stagnation  pressur  s  P,  ,  P,  ,  and  P2  are  equal 
and  can  be  denoted  as  P.  P  j 

Furthermore,  since  the  contour  of  the  mixing  chamber  in  this 
case  varies  such  as  to  maintain  the  constant  pressure 
requirement,  the  area  at  the  exit  of  the  mixing  chamber  A2 
cannot  be  equated  to  A,  t-A|s  >  but  roust  be  determined  from 
the  analysis.  ^ 

It  should  be  noted  that  if  the  flow  from  the  mixing  chamber 
is  directly  discharged  into  the  atmosphere,  the  pressure 
throughout  the  ejector  would  be  atmospheric.  In  such  a  case, 
there  would  be  no  secondary  flow  entrainment  and,  therefore, 
no  thrust  augmentation,  i.e.,  cfi  -  1.0.  It  therefore  follows 
that  a  diffuser  is  required  to  create  a  pressure  drop  in  the 
mixing  chamber  which  would  entrain  secondary  flow  in  order 
tc  produce  thrust  augmentation  (<£>I.O). 

The  idealized  analysis  for  constant  pressure  mixing  can  be 
performed  in  a  manner  similar  to  that  for  constant  area 
mixing.  Thus,  using  equation  (19),  the  relationsnip  for 
thrust  augmentation  ratio  <£  can  be  rewritten  in  the  following 
form: 
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In  order  to  obtain  the  solution  for  the  thrust  augmentation 
ratio  for  constant  pressure  mixing,  it  is  now  necessary  to 
determine  the  area  ratio  A2/Ais  as  well  as  the  velocity 
ratios  and  V,  /V|  .  This  can  be  accomplished  as 

follows:  P  P 

Using  the  continuity  equation  between  stations  (1)  and  (2), 
there  results 


AlpV,p+AlsV,s  =  A2V2 


(38) 


Also,  applying  momentum  considerations  between  stations  (1) 
and  (2)  and  considering  the  constant  pressure  condition  at 
the  mixing  chamber  walls,  there  follows 


r[z)  pAz  V2  /5A|pV,p  pA|  V| 
PAlp+PA,s-PA2-jPdAn=— - - ^ - p  09) 


where  dAn  is  the  axial  component  of  the  incremental  ejector 
surface  area. 


r 

The  integral  J  PdAn  is  given  by 
JU) 

J  PdAn=P(Alp+Ais-A2) 


(40) 


Substituting  equation  (40)  into  (39),  there  results 


A'pVip2+A'SV>S2=A2V22 


(41) 
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Squaring  equation  (38)  and  dividing  by  equation  (41)  yields 


A2  = 


lV,p+A,sVlsJ 

V'p+A«sV's 


Equation  (42)  can  be  transformed  as  follows 


(42) 


(43) 


Using  equation  (20),  the  velocity  ratio  V)s / V,  can  be 

expressed  as  follow'*  ^ 


(44) 


Substituting  equation  (43)  into  equation  (44)  and  solving  for 
V 3 / V, p  ,  there  follows 


(45) 


Also,  from  equations  (23)  and  (27)  and  the  constant  pressure 
condition,  there  results 


(46) 
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Equating  equations  (45)  and  (46)  and  simplifying,  the 
following  quadratic  equation  for  Vis/V|p  is  obtained: 


(S)\  qp  (V)  _  Vqd2-1  . 0 

V,P  aE(a0~/yao~\ )  lp 

Solving  equation  (47)  for  V)s / V( p  ,  there  results 


(47) 


(48) 


It  can  be  noted  that  the  velocity  ratio  V|S/V|p  must  be 
greater  than  zero.  Since  in  equation  (48)  the  term 


'^//otD+4aE ‘s/ qp~l  (etc-  Vop-!)  >flp 


for  aD  >  1.0  and  aE>  0,  the  velocity  ratio  V|s/Vtp  will  be 
greater  than  zero  only  if  the  positive  sign  in  front  of  the 
square  root  term  is  used.  Hence, 


Finally,  substituting  equations  (43)  and  (46)  into  equation 
(37),  the  thrust  augmentation  ratio  for  constant  pressure 
mixing  can  be  expressed  as  a  function  of  area  ratios  aE 
and  aD  and  the  velocity  ratio  V,s/V|p  .  Thus, 


(50) 
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where  Vis/V|p  ,  as  given  in  equation  (49),  is  also  a  function 
of  the  area  ratios  aE  and  ao  only. 

Similarly,  the  mass  entrainment  ratio  w  can  be  expressed  as 


/  V i s \  _  °p2~*~ ^qE'^o~^^ap~ (51) 

'V,P'  2(a0-y^7) 

For  the  special  case  of  constant  pressure  mix'  ,  ejector 
without  diffuser,  i.e.,  a0  -  1.0,  it  is  evident  from 
equation  (49)  that  V,  / Vj  =  0. 

The  above  result  indicates  that  for  the  case  of  no  diffuser, 
i.e.,  a0  1.0,  there  will  be  no  secondary  flow  and, 
therefore,  no  thrust  augmentation. 

B.  PRACTICAL  ANALYSIS 

Presented  in  this  section  is  a  practical  analysis  which 
includes  the  effects  of  many  important  jet  ejector  flow 
parameters  neglected  in  the  idealized  approach  yet  yields 
solutions  with  a  relatively  reasonable  computational  effort. 
This  analysis  pertains  to  a  jet  ejector  with  constant  area 
mixing  having  a  conical  diffuser  attached  to  the  exit  of  the 
mixing  chamber. 

As  mentioned  previously,  the  practical  analysis  is  performed 
in  the  following  stages: 

1.  The  Incompressible  Analysis  With  Flow  Losses 

Although  this  analysis  is  based  on  the  .incompressible  fluid 
flow  conditions,  it  includes  the  following  practical 
operating  conditions  and  flow  Losses: 

(a)  Effect  of  forward  speed  (parallel  to  ejector 
walls) . 

(b)  Effect  of  diffuser. 
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(c)  The  total  head  loss  at  the  secondary  entrance  to 
the  mixing  chamber. 

(d)  The  friction  head  loss  at  the  mixing  chamber  walls. 

(e)  Total  head  loss  in  the  diffuser. 

For  the  incompressible  analysis  including  the  flow  losses 
and  the  operating  conditions  specified  above,  the  jet  ejector 
flow  equations  (1)  through  (10)  can  be  reduced  to  the 
following : 

Continuity  in  the  Mixing  Chamber 


V|p+  aE  Vls  =  (aE+l)V2 


(52) 


Continuity  in  the  Diffuser 


V2  =  a0V3 


(53) 


Momentum  Across  the  Mixing  Chamber 


(aE+l)(P2-P,) 


P\2 
q  q 


(aE  +  l)pV22  (aE4l)f(l/D)p(Vls  +  V2)2 
q  2g 


(54) 


Bernoulli^  Equation  for  the  Secondary  Flow  up  to 
the  Entrance  to  the  Mixing  Chamber 


pv  oc‘ 
2q 


--P,  + 


p{!+XE)V,s 

2g 


(55) 
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* 


Bernoulli's  Equation  for  the  Flow  in  the  Diffuser 


P2  + 


1  [(aD  I )  X oIQd  I)  ]  Vj 

2i 


(56) 


The  solution  of  the  system  of  five  equations,  (52)  through 
(56),  is  similar  to  that  presented  above  for  the  idealized 
theoretical  analysis.  In  the  final  analysis,  the  five 
equations  can  be  reduced  to  one  quadratic  equation  for  the 
nondimensional  mixed  flow  velocity  at  the  ex;t  of  the 
diffuser  V3/V|p.  This  equation  is 


2  r  2  2  2' 

)+/jl  aE  (aD  — I)  +  2aEaD 


+ 


XE(aE  +  l)2a02+  f  (L/DH2aE+l)2aDZ+  X0aE2(aD-l)' 


2  2.  2.  .21 


/[f(L/DH2aE+l)V+  X0aE(a0S-l)]}  ^ 


2  2 
{2(aE-0ao-4(a  aEaQ-2[f  (l/DH2aE-H)  a0+ XE(aE-H)a0  + 


2^2f(L/D)(2aE  +  i)aD}(— - 

\vtp  / 


{ ( 2 aE  —  I )  t^2aE(aE-2)-[xE+f(L/0)l  +  ^2f(L/D)}  =  0  (57) 
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Also,  the  nondimensionalized  secondary  velocity  ratio 
v's/v'P  is  §iven  by 


(aE-H)a0^V|3J- 1 

UP  / 

aE 

(58) 


Ihe  thrust  augmentation  ratio  <£  for  the  ejector  with  forward 
speed  Vec  is  given  by 


(aE  +  l)aDV3(V3-V0C)+VlnV( 


p-a 


ap 


(59) 


Equation  (59)  can  be  expressed  in  the  following  form: 


4>  z  2 

atf  ~  (I  +  \E) 


r — ' - N7V  r2-{aeWiHl-//)aD(r^)  - 

[laE  +  l)a0(~-j-l]  ^V|P/ 


h  aE 


where 


(61) 


The  mass  entrainment  ratio  w  is 


W  =  aE 


(62) 


The  above  analysis  is  utilized  to  determine  the  effects  of 
various  flow  losses,  diffuser  area  ratios,  and  forward  speed. 
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2 .  Effect  of  Flow  Compressibility 

As  discussed  previously,  the  effect  of  flow  compressibility 
is  treated  assuming  no  losses,  no  uiffuser,  and  no  forward 
speed.  This  analysis  utilizes  the  flow  equations  (1) 
through  (10),  with  the  following  conditions: 


Vcc  =  f(L/D)  =  XE  ;XD=0 
ao  zVdz  Id  :  ‘  0 
P2  1  PQ 


(63) 


I'sing  the  above  stipulated  conditions,  the  jet  ejector  flow 
equations  (1)  through  (10)  can  now  be  reduced  to  seven 
equations  with  seven  unknowns,  which  are  /5ip.  y0|S ,  pz  ,  Vtp, 

V,s  ,  V 2  ,  and  P,  .  These  equations  in  their  non-  r 

dimensional ized  form  are  as  follows: 


Continuity  in  the  Mixing  Chamber 


Pat 


V'p\  (*>'  swV 

+  aE  y 


ap 


Js\  .  ‘  ‘  /  v2 

-  (aE  -H ) 


Pan/  \Vqp )  '  '  \Pap/  VVqp 


(64) 


Moment urn  Across  the  Mixing  Chamber 


IdfHHpQ-P,) 


P'p\  aC  /  P*s\  (  Pl\  /_v  2  2 

J  a  Ip.in/  WftJ  q  l<5fin/\v 


PapvQp  3  \p Op/  "/ap/  9  Vap/  9  'Pop'  ^v0p 

...(65) 


Conservation  of  Energy  in  the  Mixing  Chamber 


P'P 


V, 


Pop/  \vop 


P'  +  °c( : 


fP's\  /V|s\  /  Ta\  ^e  +  i)/P2\/v2 


O  \VQn  /  \POp/  V  Qp  /  VOp 


'  1  a\  vat 
Ton  /  C0T( 


-W- 


P’0pJ  \pQp/  Vvop  !  V/'1 


P2  vA 
77  +  2T) 


..(66) 
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Conservation  of  Energy  for  the  Primary  Flow  up  to 
the  Nozzle  Exit 


2gcpT0pJ  2gcpT,p  J 


+ 


Vq, 


Vfl. 


V0 


(67) 


Conservation  of  Energy  for  the  Secondary  Flow  up  to 
the  Entrance  to  the  Mixing  Chamber 


2gCpTQJ  2  qcp  T1$J 


hs 


vap  %  v  vtp ' 

Isentropic  Process  of  the  Primary  Flows 


(68) 


p>  s  r_r _ !i_i 

(y-i  CnTn.JJ 


Y-' 


Op  * x 


Pip  cPTOp^ 

Isentropic  Process  o£  the  Secondary  Flow 


P, 


J-I - — - 1 

‘•X-1  PtjCp^Q'J-1 


X 

r-< 


(69) 


(70) 


In  these  equations,  ^ap  and  Vop  refer  to  the  density  and 
velocity,  respectively,  for  the  primary  jet  exhausted 
directly  to  the  ambient  and  are  given  by 


Pac 


CP  Top 


x-i 
P\  X 


J  VP 


(71) 


x-t 


2’CPTOp4-(i£)  "  ] 


(7  2) 


31 


The  solution  of  the  system  of  equations  (64)  through  (70)  is 
performed  as  follows: 

From  equations  (69)  a^d  (70), 


(73) 


Substituting  equate  :  (72)  into  equation  (67),  there  follows 


(74) 


Since 


and 


(75) 
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Substituting  equation  (75)  into  equation  (74)  _ ields 


(76) 


Similarly,  substituting  equation  (72)  into  equation  (68), 
there  follows 


Since 


and 


(77) 


P. 


P 


(78) 
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Substituting  equation  (78)  into  equation  (77)  yeilds 


Using  the  basic  iser.-ropic  relationship, 

/  -P|p\ .  /  T'p\rrr 
\PapJ  VQp  / 

Substituting  equation  (80)  into  (73)  yields 

t  _L  o  Izi 

/T|pWp°  \  y 

/  Pis ^  VQp/  vop/ 

u"‘  ft) 


(79) 


(80) 


(81) 


From  the  equation  of  state, 


(82) 
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Also,  from  the  isentropic  relationship, 


P. 


Pa 


y 

T|p\xrr 

Tap/ 


(83) 


Substituting  equation  (82)  into  equation  (83),  there  follows 


P.  = 


rl 

y 


(84) 


Making  use  of  equation  (64)  and  substituting  equations  (76), 
(77),  (80),  (81),  (82),  and  (84)  into  equation  (65),  there 
results 


(85) 
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Similarly,  equation  (66)  can  be  transformed  as  follows: 


(86) 


Equations  (85)  and  (86)  are  the  two  final  equations  which 
contain  the  unknowns  T,p/VQp  and  Vg/Vpp  and  can  be  solved 
by  an  iteration  method.  It  is  seen  that  Cp,  the  specific 
heat  at  constant  pressure,  no  longer  appears  in  the 
equations . 

The  thermodynamic  properties  of  the  flows  at  different 
stations  along  the  ejector  can  be  expressed  in  terms  of 
T'p/Tap  anc*  V2/Vop  as  follows: 
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(95) 


(96) 


(97) 


Finally,  the  thrust  augmentation  ratio  <f>  and  the  mass 
entrainment  ratio  w  can  be  determined  as  follows: 

(98) 


(99) 

3 .  Nonuniform  Velocity  Profile  at  the  Secondary  Entrance 

In  both  the  idealized  and  practical  analyses,  a  one¬ 
dimensional  approach  is  used  which  assumes  the  secondary  flow 
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velocity  at  the  entrance  plane  as  uniform,  irrespective  of 
the  secondary-to-primary  area  ratio.  For  the  case  of  no 
diffuser,  the  idealized  analysis  yields  a  maximum  thrust 
augmentation  ratio  of  2.0  at  an  infinite  area  ratio.  This 
is  obviously  not  correct,  because  when  the  area  ratio  tends 
to  infinity,  the  ejector  configuration  reduces  to  a  free  jet, 
in  which  case  the  augmentation  ratio  must  be  unity.  It  is 
believed  that  this  discrepancy  is  caused  by  the  assumption  of 
the  uniform  flow  velocity  at  the  secondary  entrance.  In 
othe-r  words,  the  applicability  of  the  one -dimensional  approach 
is  limited  to  not  too  large  secondary-to-primary  area  ratios, 
and  a  three-dimensional  analysis  is  required  for  relatively 
large  area  ratios.  However,  such  a  task  has  never  been 
undertaken  by  any  of  the  investigators,  and  it  is  beyond  the 
scope  of  this  program  to  formulate  a  solution  to  this 
complex  problem. 

Instead,  an  empirical  correction  factor  is  herein  developed 
to  account  for  the  effect  of  the  nonuniform  secondary  flow 
velocity  profile.  The  derivation  of  this  factor  is  as 
follows : 

Under  the  assumption  of  uniform  secondary  flow  velocity,  the 
momentum  equation  for  the  control  volume  inside  the  mixing 
chamber  with  no  diffuser  and  no  losses  becomes: 


(aE+l)  (Pa-P, ) 


P 

9 


+  aEV2 


's 


(100) 


By  using  Bernoulli's  equation  for  the  secondary  flow  aud  the 
ambient  conditions  in  front  of  the  ejector,  equation  (100)  is 
reduced  to 


(101) 


In  these  equations,  the  ’’bar’'  denotes  uniform  values. 
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If  the  secondary  velocity  is  not  uniform,  the  momentum 
equation  becomes 


(aE  +  1)  Pa-P, 


P_ 

Q 


-laE  +  l)V 


2 

2 


(10  ) 


The  solution  of  equation  (102)  for  (Vg/Vip)  is  bas^d  on  the 
following  assumptions: 

(a)  The  primary  velocity  at  the  entrance  plane  is 
still  the  same  as  in  the  case  of  uniform  secondary 
velocity. 

(b)  The  primary  pressure  at  the  entrance  plane  is 
the  same  as  in  the  case  of  uniform  secondary 
velocity.  Hence,  P,p  -  P( 

(c)  The  exit  velocity  of  the  mixed  flow  V2 ,  though 
different  from  that  in  the  case  of  uniform 
secondary  velocity,  is  uniform  at  the  mixing 
chamber  e.xit. 

Since  the  secondary  flow  has  a  uniform  stagnation  pressure 
PQ  far  upstream,  by  assuming  parallel  streamlines  at  the 
entrance  plane  it  is  possible  to  evaluate  the  static  pressure 
distribution  Pi5  by  making  use  of  Bernoulli's  equation. 

Thus,  equation  (102)  is  finally  reduced  to 
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Equations  (101)  and  (103)  yield  the  exit-to-primary  velocity 
ratios  based  on  uniform  and  nonuniform  secondary  entrance 
velocity  profiles,  respectively.  The  corresponding  thrust 
augmentation  ratios  for  the  two  cases  are  given  below. 

The  thrust  augmentation  ratio  for  the  case  of  uniform 
secondary  velocity  is 


<t>  - 


(104) 


where  Va p/\  0  and  V2/V|p  can  be  obtained  from  the  idealized 
theoretical  analysis  as  follows: 


2 


Qe-I 


(105) 


and 


/VA2  f  -(aE-l)  +  aE  " 

Wip/  L  aEZ+l 


(106) 
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For  the  nonuniform  case,  the  thrust  augmentation  ratio  is 


4>* 


(107) 


_  2 

where  (V2 / V, p )  as  given  by  equation  (103)  will  be 
determined  below. 

The  correction  factor  for  thrust  augmentation  ratio  due  to 
nonuniform  velocity  profile  at  the  secondary  entrance  is 
herein  defined  as 


<*>  W.p  / 

f '  /  vj_y 

\V,p  / 


(108) 


It  is  now  necessary  to  establish  relationships  for  V,s/V|p 
and  V,s/V,p  to  determine  V2/V,p  given  by  equation  (103). 
According  to  Reference  26,  the  velocity  profile  of  a  flow 
entrained  by  a  free  jet  can  be  well  approximated  by  a 
cosine  curve.  Therefore,  it  is  reasonable  to  assume  that  a 
cosine  curve  distribution,  as  shown  in  Figure  2,  will 
adequately  represent  the  velocity  profile  at  the  secondary 
entrance  of  a  jet  ejector  system. 

It  is  further  assumed  that  the  local  secondary  velocity  at 
the  Jet  periphery  (a)  is  the  same  as  the  uniform  velocity 
V|S  .  Thus,  the  secondary  entrance  velocity  profile  at  any 
radial  station  can  be  expressed  as  follows: 


I  +  cos 


(r‘-l)  w 

K 


z 


(109) 
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In  equation  (109) „  r'  is  the  radial  distance  from  the  ejector 
center  line,  nondimensionalized  by  the  jet  radius,  and  xy  is 
a  constant  defining  flow  nonuniformity  at  the  secondary 
entrance.  This  equation  applies  only  for  the  range 
0  <  ( r I )  / <  1.  In  other  words,  as  shown  in  Figure  2,  the 

local  velocity  V1$  gradually  decreases  to  zero  at  point  (b) 
and  maintains  zero  at  any  larger  radial  distance  from  the 
ejector  center  line.  The  location  of  oint(b) depends  on  the 
constant  k. 

By  assigning  a  suitable  value  of  k ,  the  integral 
/  (V,s/V|p)2  dA|S/A|S  in  equation  (103)  can  be  evaluated  for  all 
values  of  aE  in  terms  of  V)s  .  Let 

f  2 

/V|cd  A|c 

g  =  L . ig _ !i_  (nc 

v'sA>s 

Substituting  equation  (109)  into  equation  (110),  there 
results 


.  aE+l  /•' f  (r‘-l)ir 
O_L_jJ1  +  c0s__ 


r'd'r 


Integrating  equation  (111)  with  respect  to  x  yields 


(111) 


e-r+ 


3  1  .  (VaE+l  -Dir 


sin 


8  4aE-?r 

2  r  (vSi+l-Dir 


4  +  cos 


I  - 


cos 


8aETr 

for  aE  <  k  (k+2)  and 


][ 


9  +  cos 


(./aE+l  -I)tt 
k 

(VSi+T-Dir] 


]- 


(112) 


Zk(k +  2)tt  -16k 
8aE7r2 


(113) 


for  Qr  <  k(k  +  2) 
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The  correction  factor  fcr  the  thrust  augmentation  ratio  given 
by  equation  (108)  thus  becomes 


X  = 


(114) 


where  £  is  given  by  equation  (112)  or  (113)  and  V|5/V|D  can 
be  obtained  from  the  idealized  analysis  equation  (32)  as 
follows : 


V|s  .  [~(Qe~0  4-  ag^/2a£]-(ag-f  I)  (115) 

Vip  aE(aE2+l) 

Equation  (114)  is  utilized  to  compute  the  thrust  augmentation 
ratio  as  affected  by  the  nonuniform  velocity  distribution 
at  the  secondary  entrance.  The  results  are  herein  presented 
in  Figure  3,  which  shows  the  variation  of  <fc  vs.  aE  for 
constant  values  of  the  parameter  k.  By  examining  this 
figure,  it  can  be  noted  that  as  k  tends  to  infinity,  the  <£ 
vs.  aE  relationship  reduces  to  that  as  obtained  by  the 
idealized  analysis  with  uniform  secondary  velocity  profile. 
Also,  for  any  finite  value  of  k,  <fc  reaches  a  maximum  value 
(less  than  2.0)  at  a  finite  qee  and  then  decreases  to  unity 
as  aE  tends  to  infinity. 

In  order  to  determine  the  parameter  k ,  appropriate  test  data 
are  required  which  could  be  used  to  determine  the  area  ratio 
at e  at  which  the  thrust  augmentation  ratio  is  maximum. 

However,  most  of  the  available  test  data  (see  Section  VI)  are 
obtained  for  ejectors  having  area  ratios  aE  of  less  than  100, 
at  which  values  the  thrust  augmentation  ratio  <f>  does  not 
reach  a  maximum  value.  For  example,  the  correlation  shown  in 
Table  II  (Section  VII)  indicates  that  at  an  area  ratio  of 
aE  -  104,  the  rate  of  increase  of  <f>  with  aE  is  quite  small. 

It  is  reasonable  to  assume,  therefore,  that  for  aE  ranging 
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from  150  to  200,  the  thrust  augmentation  would  have  reached 
its  peak  value.  Therefore,  it  follows  from  Figure  3  that  a 
practical  value  of  the  parameter  k  would  b°  between  k  -  60 
and  k  -  70.  This,  however,  necessitates  further  experi¬ 
mental  evidence. 

4.  Analysis  of  Various  Elector  Types 

The  most  common  types  of  jet  ejectors  consist  of  either 
single  nozzle,  multiple  nozzle,  or  annular  nozzle  configu¬ 
rations.  These  configurations  differ  from  each  other 
primarily  in  the  efficiency  at  which  the  primary,  low 
momentum  flow  is  converted  into  a  high  momentum  flow.  The 
conversion  takes  place  in  the  mixing  chamber,  which,  in  order 
to  reduce  wall  friction  losses,  should  be  of  minimum  length. 
No  precise  method  exists  at  present  to  determine  the  minimum 
mixing  chamber  length  required  to  achieve  complete  mixing 
between  the  primary  and  secondary  flows.  The  following 
analysis  will  provide,  however,  a  first-order  approximation 
of  the  relative  effects  of  the  different  configurations  on 
the  mixing  chamber  length. 

First,  it  is  assumed  that  the  mixi  g  chamber  consists  of  two 
parts  as  shown  in  Figure  4.  The  first  part  (Lt)  ,  which 
allows  for  the  primary  jet  expansion,  starts  at  the  jet  exit 
and  ends  at  a  point  where  the  jet  boundary  reaches  the 
ejector  walls.  The  second  part  (L2)  starts  from  thereon  and 
extends  to  the  exit  of  the  mixing  chamber,  where  the  mixing 
process  is  assumed  to  be  completed. 

a.  Single  Nozzle  Configuration 

The  work  of  Squire  and  Trouncer  (Reference  19)  is 
utilized  to  determine  the  length  (Lt)  for  a  single¬ 
nozzle  primary  flow  configuration. 

In  this  reference  an  analysis  is  performed  for 
determining  the  expansion  angle  of  a  single, 
circular  jet  discharging  into  a  free  stream  of 
uniform  velocity.  If  it  is  assumed  that  this 
expansion  angle  is  not  affected  by  the  presence  of 
the  mixing  chamber  walls,  the  primary  jet  boundaries 
can  be  plotted  as  a  function  of  the  secondary  flow 
velocity  as  shown  in  Figure  5.  In  this  figure  the 
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(c) 


f - 

(d) 


Figure  4.  Definition  of  Paraneters  for 
Single  rnd  Annular  Nozsle  Configurations. 
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Figure  5.  Jet  Boundaries 


strearowise  distance  from  the  primary  jet  exit  is 
denoted  by  X,  and  the  radial  distance  of  the  jet 
boundary  from  the  jet  center  line  is  denoted  by  r. 
Typical  variations  of  V|S/Vtp  versus  secondary-to- 
primary  area  ratio  aE  are  snown  in  Figure  6. 

Using  Figures  5  and  6,  the  length  L|  can  now  be 
determined  as  follows: 

(i)  A  value  aE  is  assumed,  and  from 

Figure  6  the  corresponding  value  of 
V|S / V|p  is  obtained. 

(ii)  The  value  of  r/a  is  obtained  from  the 
following  relationship: 


r/a  =7^7+7  (U6) 


(iii)  With  the  above  values  of  V|E/Vjp  and  r/a, 
the  corresponding  value  of  X/a  is 
obtained  from  Figure  5. 

(iv)  Finally, 


±±:±(X/£\  (117) 

0  2  \  r  /a  / 

The  above  calculations  were  repeated  for  several 
values  of  aE ,  and  the  results  are  plotted  in 
Figure  7  for  the  incompressible  and  the  compressible 
analyses.  Similar  curves  can  also  be  obtained 
including  various  flow  losses.  Figure  7  indicates 
that,  for  practical  values  of  aE ,  the  ratio  (L,/D) 
is  of  the  order  of  3  to  4. 

Insufficient  data  exist  to  determine  the  length  of 
the  second  part  of  the  mixing  chamber,  L2  • 

However,  it  will  be  assumed  that  this  length  is 
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Figure  7.  Variation  of  Length  of  the  First  Part 
of  Mixing  Chamber  with  Secondary- to- Primary 
Area  Ratio  (Ho  Diffuser;  Mo  Forward  Speed; 

Do  Losses). 
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proportional  to  the  contact  surface  area  between  the 
primary  and  the  surrounding  flows. 

b.  Multiple  Nozzle  Configuration 

The  first  portion  of  the  mixing  chamber  length,  L,  , 
for  a  multiple  nozzle  configuration  can  easily  be 
determined  for  the  two-dimensional  ejector  by 
assuming  that  these  nozzles  are  evenly  spaced  with 
respect  to  the  mixing  chamber  and  with  respect  to 
each  other. 

For  a  multiple  nozzle  configuration  having  N  number 
of  nozzles  of  the  same  total  area  as  the  single 
central  nozzle,  (Aip)s  ,  the  primary  area  of  each 
nozzle  is 


N 


(A,J 


p's 


(118) 


Also,  if  the  total  secondary  area  of  the  multiple 
and  single  nozzle  configurations  is  the  same,  it 
follows  that 


(A,sV  N  (Vs 


(119) 


Therefore,  from  equations  (118)  and  (119)  it  follows 
that  the  secondary-to-primary  area  ratio  for  each 
nozzle  of  multiple  nozzle  configuration  is  the  same 
as  that  for  the  corresponding  single  central  nozzle 
type,  i.e., 


(a E)  : 

Em  (V- 


(A,J, 


(A^) 


p's 


(120) 
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Since  the  area  ratios  are  the  same  (equation  120), 
using  Figure  7,  the  effective  (L t /D )  for  each 
nozzle  of  the  multiple  nozzle  configuration  is  the 
same  as  that  for  the  single  nozzle  design,  i.e., 


(121) 


From  equation  (119),  the  diameter  of  the  effective 
mixing  chamber  for  each  nozzle  of  the  multiple 
nozzle  configuration  can  be  determined  as  follows: 


2  2 
7tD  \  _  J_  / 7rP\ 

k  4  /M  N  \  4  /S 


(122) 


or 


i(V  7f 1D’5 


(123) 


From  equations  (121)  and  (123),  it  follows  that  the 
mixing  chamber  length  (L|)  required  for  the  jet 
boundary  to  reach  the  ejector  walls  for  the  multiple 
nozzle  configuration  is  reduced  by  a  factor  of 
iA/N~  of  the  corresponding  single  nozzle  length. 
Thus , 


<12< 

Assuming  that  thf  second  part  of  the  mixing  chamber 
length  L2  is  a  ^.unction  of  the  total  contact  area 
of  the  primary  flow  with  the  surrounding  air,  the 
re1 fltionship  for  the  mixing  chamber  lengths  (L2)m 
and  _2)s  of  multiple  and  single  nozzle  configura¬ 
tions  can  be  obtained  as  follows: 
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The  contact  area  of  the  prir^ry  jet  with  the 
surrounding  flew  for  a  single  central  nozzle  type 
is  given  by 


(S)s-  7t(D)s(L2)s 


(125) 


Similarly,  the  total  contact  area  of  the  primary 
jet  for  a  multiple  nozzle  configuration  of  N 
nozzles  is 


( S )M  =  N  7T  (D)m(L2)m 


(126) 


Substituting  equation  (123)  into  (126)  yields 


(S)m  =  7N  tt(D)s(L2)m 


(127) 


Assuming  that  in  each  case  the  contact  areas  are 
such  as  to  produce  complete  mixing  at  the  exit  of 
the  mixing  chamber,  equations  (125)  and  (127)  yield 


(128) 

Combining  equations  (124)  and  (123),  there  follows 


1lw7n  (l)s  (129) 
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Equations  (124),  (128),  and  (129)  indicate  that  for 
a  multiple  nozzle  jet  ejector  configuration,  each 
part  of  the  mixing  chamber  length  as  well  as  the 
total  length  reduces  by  a  factor  of  IA/InT  as 
compared  with  the  single  central  nozzle  configura¬ 
tion.  This  reduction  in  the  mixing  chamber  length 
results  in  the  reduction  of  head  loss  due  to  wall 
frictio. .  Therefore,  the  multiple  nozzle  configu¬ 
ration  would  yield  a  superior  performance  as 
compared  to  that  of  the  equivalent  single  nozzle 
ejector  if  the  complete  mixing  conditions  at  the 
exit  of  the  mixing  chamber  are  satisfied  in  each 
case . 

c .  Annular  Nozzle  Configuration 

For  the  case  of  an  annular  primary  jet  located 
midway  between  the  walls  and  the  center  of  the 
mixing  chamber  walls,  both  parts  of  the  mixing 
chamber  lengths  (Lt)A  and  (L2)a  >  as  discussed  above, 
will  also  substantially  reduce  as  compared  to  the 
equivalent  single  nozzle  configuration. 

The  analysis  for  determining  the  mixing  chamber 
lengths  for  the  annular  nozzle  ejector  can  be 
performed  as  follows: 

Using  Figure  4(a;,  the  primary  area  of  the  single 
nozzle  configuration  is  given  by 


(A'p)s=fd' 


(130) 


From  Figure  4(b),  the  corresponding  primary  area 
for  the  annular  nozzle  is 


(SVTUe-d.2) 


(131) 
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For  Che  same  primary  jet  areas,  equations  (130) 
and  (131)  yield 


d2  =  d| -df  =  (d2  +  ti, )  (d2-d,) 


(132) 


Since  the  annular  jet  is  centrally  located  within 
the  mixing  chamber  having  a  fixed  diameter  D,  it 
follows  that 


(133) 

2  2 


Equations  (132)  and  (133)  yield 


•j 


2 


(134) 


Since  the  primary  jet  areas  for  the  two  configura¬ 
tions  are  the  same  and  since  outside  diameters  D 
of  the  mixing  chambers  are  also  the  same,  it 
follows  that 


lVs’(VA  =  T(DZ-d')  (135) 

From  the  above  analysis  it  therefore  follows,  that 
the  secondary-to-primary  area  ratio  of  the  annular 
nozzle  ejector  is  equal  to  that  of  the  equivalent 
single  central  nozzle  configuration,  i.e., 


Using  equations  (130)  and  (135),  the  corresponding 
secondary-to-primary  area  ratios  are 


aE  = 


Ais  D2 

Aip 


(136) 


(137) 


Assuming  identical  primary  jet  stagnation  conditions, 
the  velocity  ratios  Vis/V  0  will  also  be  equal. 

Using  Figure  5,  it  therefore  follows  that  the 
primary  jet  expansion  angle, 


6  -  tan 


s 


« 


for  annular  nozzle  configuration  and  the 
corresponding  single  jet  will  be  the  same.  Thus, 
using  Figure  4(c),  the  first  part  of  the  mixing 
chamber  length  for  a  single  nozzle  configuration 
can  be  expressed  as  follows: 


(L  i)$  = 


D-d 

2 


Ton d 


(138) 


Similarly,  using  Figure  4(d),  the  corresponding  J 

length  for  the  annular  nozzle  ejector  is 


(L.)a= 


0-(d2-d,) 

4 


tan  6 


(139) 


1 
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Equations  (138)  and  (139)  yield 


5,)a  _  i  r  D-tdjj-d,)-] 

(L,)$  '  2  L  D-d  . 


(140) 


Substituting  equation  (134)  into  equation  (140), 
there  follows 


Ma  I  d 

"(L|)s  1  2  li+  D  } 


(141) 


Substituting  equation  (137)  into  (141)  yields 


(L,)s  2  \  v^+T / 


(142) 


or 


<L,,»’t(l+'^f),L^  (143: 

Equation  (143)  shows  that  the  first  part  of  the 

mixing  length  (Li  )A  for  the  annular  nozzle 

configuration  is  a  function  of  the  area  ratio  aE 

and  that  for  large  values  of  aE  it  is  approximately 

equal  to  a  half  that  of  the  corresponding  single 

central  nozzle  ejector.  For  small  values  of  the 

area  ratios,  this  length  approaches  ( L , )  . 

s 

The  second  part  of  the  mixing  chamber  length  for  the 
annular  nozzle  ejector  can  be  determined  as  follows: 


The  surface  contact  area  of  the  primary  jet  with 
the  surrounding  flow  for  a  single  nozzle  ejecto*.  is 


(S)s=  7rd(L2)s  (144) 


The  corresponding  contact  area  for  the  annular  jet 
is  given  by 


(S)A=  7r(d2  +  d,)(L2)A 


(145) 


Assuming  that  in  each  case  the  seme  constant  contact 
surface  area  is  required  to  complete  the  mixing 
process  at  the  exit  of  the  mixing  chamber, 
equations  (144)  and  (145)  yield 


s 


d2+d, 


(L»), 


2'S 


(146) 


Substituting  equation  (1*3)  into  (146)  yields 


(l2]a=  'd”{l2) 


(147) 


From  equations  (142)  and  (147),,  tu.cre  results 


I 

VaE4- 1 


(148) 


Equation  (148)  shows  that  the  second  part  of  the 
mixing  chamber  length  for  the  annular  nozzle 
configuration  reduces  by  a  factor  of  I /VaE  + 1 
as  compared  to  that  of  the  equivalent  single  jet 
ejector. 
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Finally,  using  equations  (143)  and  (148),  the 
total  mixing  length  for  the  annular  jet  ejector 
is  given  by 


V.  EFFECT  OF  PARAMETERS  ON  JET  EJECTOR  PERFORMANCE 

Presented  in  this  section  are  the  results  of  a  study  showing 
the  effects  of  various  aerodynamic,  thermodynamic,  and 
geometric  flow  parameters  on  ejector  thrust  augmentation 
performance.  The  study  includes  the  following  parameters: 

A.  FLOW  LOSSES 

The  effect  of  flow  losses  is  to  reduce  ejector  performance. 
This  effect  can  oe  < learly  seen  frcm  the  nomographs  of 
Section  VII  and  requires  no  further  explanation. 

B.  DIFFUSER 

The  results  of  the  idealized  analysis  presented  in  Figures  8 
and  9  show  that  under  static  conditions,  at  any  secondary- to¬ 
primary  area  ratio,  an  increase  in  diffuser  area  ratio  yields 
better  thrust  augmentation  performance.  However,  as  shown  in 
Figure  10,  this  is  not  true  for  a  practical  case  when  the 
flow  losses  are  included.  Although  the  overall  diffuser  loss 
factor  XD  varies  with  the  area  ratio  aD,  for  the  purpose  nf 
this  duscussion  all  flow  losses  including  the  diffuser  loss 
factor  XD  are  held  constant  (at  typical  values)  for  all 
diffuser  area  ratios. 

Examining  Figure  10,  it  can  be  seen  that  for  any  pracitcal 
value  of  secondary-to-prima.y  area  ratio  aE ,  the  thrust 
augmentation  ratio  reaches  an  optimum  value  with  the 
diffuser  area  ratio  a0  ranging  between  1.5  and  2.0.  For  the 
diffuser  area  ratios  larger  than  2.0,  the  thrust  augmentation 
ratio  would  be  actually  lower  than  that  indicated  in 
Figure  10  as  a  result  of  increase  in  the  diffuser  loss  factor 
X o  with  increase  of  ay 

C.  FORWARD  SPEED 

Figure  11  shows  the  effect  of  forward  speed  on  thrust 
augmentation  ratio  for  an  idealized  condition  with  no  flow 
losses  and  no  diffuser.  It  can  be  seen  from  this  figure 
that  an  increase  in  forward  speed  results  in  a  substantial 
reduction  of  thrust  augr  jntation  ratio  for  all  seconder y- to¬ 
primary  area  ratios, a£.  Furthermore,  for  any  specific 
nonzero  value  of  forward  speed,  the  thrust  augmentation 
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Figure  10.  Effect  of  Diffuser  on  Thrust  Augmenta¬ 
tion  of  «  Jet  Ejector  with  Constant  Area  Mixing 
Chamber  (Ho  Forward  Speed,  With  Typical  Flow 
Losses). 
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ratio  reaches  an  optimum  value  cf  1.0  <  <  2.0,  and  then 

reduces  to  unity  as  the  area  ratio  further  increases  to 
infinity.  From  a  physical  point  of  view,  the  case  of 
infinite  aE  corresponds  to  the  condition  of  primary  flow 
discharged  directly  into  the  ambient,  resulting  in  thrust 
augmentation  ratio  of  unity.  This  figure  also  shows  the 
locus  of  maximum  thrust  augmentation  ratio  (dotted  line)  as 
a  function  of  aE  for  the  range  of  forward  speed  considered. 

Figure  12  shows  the  effect  of  forward  speed  on  thrust 
augmentation  ratio  including  typical  flow  losses  for  the 
case  of  ajj  --  2.0.  Comparing  Figures  11  and  12,  it  can  be 
noted  that  the  effect  of  forward  speed  on  the  reduction  of 
thrust  augmentation  ratio  is  more  severe  when  the  major  flow 
losses  are  included.  Specifically,  fcr  forward  speeds 
greater  than  fi  =  0.05,  the  thrust  augmentation  ratio 
continuously  reduces  with  an  increase  in  aE  and  eventually 
reaches  a  value  of  less  than  1.0  in  the  practical  design 
range  of  aE. 


The  above  discussion  refers  to  an  ejector  forward  speed  along 
the  ejector  longitudinal  axis,  i.e.,  parallel  to  jet  thrust. 
If,  however,  the  forward  motion  is  normal  to  the  ejector 
axis,  under  idealized  flow  condition,  the  ejector  performance 
should  not  be  affected  by  such  motion.  In  a  practical  case, 
the  forward  motion  of  an  ejector  perpendicular  to  the 
direction  of  thrust  will  result  in  an  asymmetric  inflow  at 
the  secondary  entrance,  thereby  causing  an  increase  in  the 
secondary  entrance  head  loss  }  __ . 

D.  FLOW  COMPRESSIBILITY 

The  effect  of  flow  compressibility  on  jet  ejector  performance 
is  herein  determined  assuming  no  forward  speed,  no  diffuser, 
and  no  flow  losses.  Typical  computer  results  for  this 
analysis  are  presented  in  Figures  13  through  15.  In  each  of 
these  figures,  the  thrust  augmentation  ratio  obtained  from 
the  idealized  incompressible  analysis  is  also  shown  for  the 
purpose  of  comparison.  It  is  seen  that  in  general  the  thrust 
augmentation  ratio  is  reduced  when  the  flow  compressibility 
is  accounted  for. 
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Ratio  (with  Diffuser;  with  Typical  Losses). 


2.0 


e  13.  Effect  <f  Primary  Stagnation  Pressure  on  Thrust 
notation  Ratio  (Mo  Forward  Speed;  Ho  Diffuser;  Mo  Losses) 
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Incompressible  Analysis 


e  IS.  Effect  of  Ambient  Temperature  on  Thrust  Augments 
Ratio  (No  Forward  Speed;  No  Diffuser;  No  Losses). 


Figure  13  shows  that  the  increase  of  the  primary  stagnation 
pressure  P0.  from  20  psia  to  25  psia  results  in  about  3 
percent  reduction  of  the  thrust  augmentation  ratio.  Pop  is 
limited  to  25  psia  in  the  calculations  in  order  to  avoid  the 
choking  condition  of  the  nozzle. 

The  effect  of  the  variation  of  the  primary  stagnation  temper¬ 
ature  'rop  is  illustrated  in  Figure  14.  The  temperature  range 
unde,  consideration  is  between  600°F  and  1400°F.  It  is  seen 
from  this  figure  that  an  increase  of  temperature  Top 
decreases  the  thrust  augmentation  ratio  up  to  a  maximum  of 
about  4  percent  at  area  ratios  around  10. 

The  increase  in  ambient  temperature,  as  shown  in  Figure  15, 
results  in  an  increase  of  thrust  augmentation  ratio  <f>  .  This 
increase,  however,  does  not  exceed  2  percent  for  the  range  of 
ambient  temperatures  raised  from  0°F  to  200°F  and  the  range 
of  <2e  considered. 

In  general,  it  can  be  concluded  that  although  thr;  effect  of 
flow  compressibility  is  to  reduce  the  jet  ejector  performance 
as  compared  to  the  idealized  flow  conditions,  the  effect  of 
the  variation  of  the  compressible  flow  parameters,  such  as 
Pop*  Top  ,  and  Tq,  seems  to  be  of  little  consequence. 

E.  MIXING  CHAMBER  SHAPE 

In  the  present  investigction,  two  mixing  chamber  shapes  have 
been  considered;  i.e.,  the  constant  area  and  the  constant 
pressure  mixing  configurations. 

The  analysis  of  the  mixing  chamber  shapes,  other  than  the 
two  investigated,  is  indeed  a  very  difficult  task  and  has 
not  been  successfully  attempted  by  any  of  the  Investigators 
in  the  past.  It  appears,  however,  that  due  to  possible  flow 
separation  at  the  mixing  chamber  walls,  a  divergent  mixing 
chamber  shape  would  not  be  suitable. 

Some  of  the  information  on  the  mixing  chamber  shapes  can  be 
obtained  by  comparing  Figures  8  and  9.  Examining  these 
figures,  it  can  be  seen  that  for  the  same  secondary-to¬ 
primary  area  ratio  and  diffuser  area  ratio  aD,  the  constant 
area  mixing  chamber  shape  yields  superior  ti  rust  augmentation 
performance  as  compared  to  that  of  a  constant  pressure 
mixing  configuration. 


F.  MIXING  CHAMBER  LENGTH 


The  data  presented  in  this  report  are  valid  only  for  the 
case  when  the  velocity  of  the  mixed  flows  at  the  exit  of  the 
mixing  chamber  is  uniform.  This  implies  complete  mixing  of 
the  primary  and  secondary  flows  which  necessitates  adequate 
mixing  chamber  length. 

The  exact  mixing  process  inside  a  mixing  chamber  of  an 
ejector  is  not  yet  clearly  understood.  Mikhail  (Reference 
20)  made  an  attempt  to  solve  the  problem;  however,  the  final 
equations  are  developed  in  terms  of  "mixing  length  parameters" 
which  vary  from  case  to  case  and  have  to  be  determined 
empirically. 

The  experimental  data  of  Reference  15  show  that  for  the 
mixing  chamber  length  of  L/D  =  12,  the  velocity  distribution 
at  the  exit  of  the  mixing  chamber  is  uniform.  The  value  of 
L/D  =  12  is  not  necessarily  representative  for  a  variety  of 
jet  ejector  configurations,  since  this  value  will  in  general 
depend  on  the  secondary-tc-primary  area  ratio  aE  and  other 
geometric  parameters. 

On  the  other  hand,  the  experimental  data  of  Reference  2 
indicate  that  an  increase  in  the  mixing  chamber  length 
beyond  the  value  of  L/D  =  5.0  or  6.0  causes  an  adverse  effect 
on  the  jet  ejector  performance.  This  implies  that  the  flow 
losses  due  to  partial  mixing,  associated  with  short  mixing 
chamber  lengths,  are  less  predominant  than  the  friction 
losses  caused  by  increasing  mixing  chamber  length  to  achieve 
complete  mixing.  It  appears,  therefore,  that  in  selecting 
practical  mixing  chamber  lengths  the  considerations  of  the 
flow  losses  due  to  partial  mixing  may  be  of  secondary 
importance. 

This  subject  is  further  discussed  later  in  the  text.  Also, 
the  charts  for  estimating  mixing  chamber  lengths  of  various 
ejector  configurations  are  presented  in  Section  VII.  These 
charts  are  based  on  the  fact  that  the  practical  mixing 
chamber  length  required  for  optimum  ejector  performance  is 
less  than  that  required  for  complete  flow  mixing.  Further¬ 
more,  it  is  assumed  that  the  flow  losses  due  to  partial 
mixing  can  be  neglected  provided  that  the  total  mixing 
chamber  length  for  a  single  nozzle  Jet  ejector  is  greater 
than  L/D  =  6.0 . 
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G.  DIVERGENCE  AND  INCLINATION  OF  NOZZLES 

In  the  previous  discussions,  it  was  assumed  that,  irrespec¬ 
tive  of  the  ejector  type,  the  orientation  of  the  primary  jet 
nozzles  was  along  the  longitudinal  axis  of  the  ejector. 

For  the  case  of  the  central  nozzle  or  multinozzles  installed 
at  an  angle  with  the  ejector  axis,  or  for  the  case  r f  the 
annular  nozzle  with  the  jet  diverging  from  the  ejector  axis, 
the  one-dimensional  analysis  is  considered  to  be  inadequate 
in  predicting  the  performance  of  such  ejector  configurations. 
The  major  effect  will  be  the  increase  of  mixing  losses  due  to 
the  flow  interaction  which  is  difficult  to  treat  analytically. 

In  the  muiticcncer.tr'  c  annular  nozzle  configuration  as 
developed  by  Bertin,  Reference  10,  the  divergence  of  the 
nozzle  seems  to  be  of  merit  insofar  as  mixing  is  concerned. 
However,  no  analytical  work,  besides  the  usual  one- dimensional 
approach,  has  been  presented  by  Bertin  for  the  purpose  of 
analyzing  the  performance  of  jet  ejectors  with  divergent 
multiconcentric  annular  nozzles. 

H.  EFFECT  OF  GROUND  PROXIMITY 

Available  test  data,  e.g..  References  21  and  22,  indicate  that 
in  the  ground  proximity,  the  thrust  augmentation  of  an  ejector 
is  reduced.  This  is  attributed  to  the  reduction  of  static 
pressure  prevailing  in  the  flow  from  the  ejector  exit.  The 
test  data  also  show  that  the  ejectors  operating  in  parallel, 
with  a  common  baseplate,  yield  a  higher  thrust  augmentation 
than  the  ejectors  operating  individually.  This  is  due  to  the 
building  up  of  pressure  underneath  the  baseplate  where  a 
stagnation  condition  exists.  This  effect  is  similar  to  what 
prevails  in  the  case  of  annular  nozzles  in  ground  proximity, 
for  which  some  theoretical  and  experimental  data  are 
presented  in  References  23  and  24. 

I .  NONUN I FORM  VELOCITY  AT  THE  SECONDARY  ENTRANCE 

The  practical  analysis  presented  in  Section  IV  as  well  as 
other  analyses  available  from  the  existing  literature,  is 
performed  with  a  one-dimensional  flow  approach.  The  necessary 
assumptions,  among  others,  are  that  the  velocities  at  the 
secondary  entrance  and  also  at  the  mixing  chamber  exit  are 
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uniform.  In  practical  operations  these  idealized  conditions 
do  not  always  apply. 

For  small  secondary- to-primary  area  ratios,  the  velocity  at 
the  secondary  entrance  is  practically  uniform.  However,  as 
the  area  ratio  aE  increases,  the  secondary  entrance  velocity 
close  to  the  mixing  chamber  walls  becomes  much  lower  than 
that  close  to  the  primary  jet. 

For  very  large  area  ratios  (as  Oe  tends  to  infinity),  the 
velocity  of  the  secondary  flow  becomes  zero  at  a  finite 
radial  distance  from  the  primary  jet.  In  fact,  this 
condition  corresponds  to  a  free  jet  discharged  directly  into 
the  ambient,  in  which  case  no  thrust  augmentation  is  possible. 
The  available  analyses,  utilizing  the  assumption  of  uniform 
secondary  flow  velocity,  result,  however,  in  a  thrust 
augmentation  ratio  of  2.0  as  aE  tends  to  infinity. 

In  the  previous  section,  an  attempt  is  made  to  account  for 
the  radial  variation  of  the  local  velocity  at  the  secondary 
entrance.  However,  this  analysis  is  based  on  a  semi- 
empirical  approach  requiring  furhter  experimental  verifica¬ 
tion. 

J.  INCOMPLETE  MIXING 

The  effect  of  incomplete  (partial)  mixing  of  the  primary  and 
secondary  flows  at  the  exit  of  the  mixing  chamber  has  also 
been  investigated.  As  discussed  previously,  the  incomplete 
mixing  will  exist  due  to  the  insufficient  length  of  the 
mixing  chamber,  thereby  resulting  in  a  nonuniform  velocity 
profile  at  the  exit  of  the  mixing  chamoer. 

The  difficulties  of  analyzing  partial  mixing  are  in  principle 
similar  to  those  associated  with  nonuniform  velocity  at  the 
secondary  entrance  discussed  above.  In  this  case,  however, 
even  if  the  velocity  profile  at  the  exit  of  the  mixing 
chamber  s  known  or  assumed,  the  corresponding  variation  of 
the  static  pressure  distribution  cannot  be  analytically 
determined. 

Forstall  and  Shapiro  (Reference  25)  indicate  that  the 
velocity  profile  at  the  exit  of  the  mixing  chamber  can  be 
quite  accurately  represented  by  a  cosine  or  three-halves 
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power  curve.  However,  information  on  the  static  pressure 
distribution  at  the  exit  of  the  mixing  chamber  and  the  mass 
flow  rate  of  the  partially  mixed  flow  is  not  available. 

Thus,  in  order  to  predict  reliably  the  effects  of  the 
incomplete  (partial)  mixing,  an  experimental  investigation  is 
necessary  to  obtain  the  pertinent  data  on  pressure  and 
velocity  distributions  at  the  exit  of  the  mixing  chamber. 

K.  MIXING  LOSSES  IN  THE  MIXING  CHAMBER 

The  mixing  losses  in  the  mixing  chamber  cannot  be  considered 
as  an  additional  external  force  which  could  be  readily 
introduced  in  the  momentum  equation.  One  method  of 
accounting  for  the  mixing  losses  is  to  assume  a  mixing 
efficiency  factor  (t]h)  in  the  flow  momentum  considerations. 
However,  there  is  no  information  available  for  determining 
this  mixing  efficiency  factor. 

L.  FLOW  VISCOSITY 

A  real  fluid  flow  is  generally  treated  as  an  inviscid  flow, 
except  that  at  the  locations  close  to  the  walls  a  boundary 
layer  effect  is  taken  into  consideration.  The  flow  viscosity 
effects  give  rise  to  friction  drag  which  in  turn  reduces  the 
jet  ejector  thrust.  In  the  present  investigation,  the  losses 
due  to  the  friction  drag  along  the  mixing  chamber  walls  are 
included  in  the  friction  factor  f,  which  appears  in  the 
momentum  equation. 

For  the  diffuser  section,  the  losses  due  to  friction  drag  and 
flow  separation  are  included  in  the  diffuser  loss  factor  X0  . 

M.  HEAT  CONDUCTION  LOSSES 

In  the  practical  analysis  (compressible  flow),  the  heat 
conduction  losses  will  affect  the  energy  equations  presented 
in  Section  IV.  However,  the  information  pertaining  to  the 
rate  of  heat  conduction  which  varies  with  ejector  geometry, 
conductivity  of  mixing  chamber  walls,  and  the  operating 
conditions  cr  the  ejector  system  is  inadequate  for  reasonable 
evaluation. 
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Furthermore,  due  to  the  heat  conduction  losses,  the  ejector 
flows  become  neither  isentropic  nor  adiabatic,  resulting  in 
an  additional  difficulty  in  analyzing  the  flow  process. 

In  general,  however,  the  mixing  chamber  is  usually  not  very 
long  and  the  secondly  flow  is  almost  under  ambient 
conditions,  so  that  heat  conduction  losses,  if  any,  through 
the  mixing  chamber  and  diffuser  walls  are  not  significant 
and  can  be  neglected  without  materially  affecting  the 
accuracy  of  the  analysis. 

N.  EFFECT  OF  SPECIFIC  HEAT  RATIO  / 

Using  th">  computer  program,  Appendix  I,  an  investigation  was 
performed  to  determine  the  effects  of  the  specific  heat 
ratio  y  on  jet  ejector  performance.  It  was  found,  however, 
that  for  a  range  of  y  from  1.25  to  1.4,  the  thrust  augmenta¬ 
tion  ratio  is  practically  unchanged  ^ss  than  0.1  percent). 

It  appears,  therefore,  that  the  effect  of  the  specific  heat 
ratio  can  be  neglected  in  the  practical  jet  ejector  analysis. 
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VI.  CORRELATION  with  experimental  data 

Available  experimental  data  which  are  suitable  for 
correlation  with  the  results  of  the  present  analysis  are 
rather  limited.  Those  data  which  contain  information  on 
thrust  augmentation  do  not,  in  general,  furnish  sufficient 
design  and  construction  details  which  could  be  used  to 
determine  the  various  flow  loss  factors.  Consequently,  for 
some  of  the  jet  ejector  configurations,  these  loss  factors 
were  determined  by  using  good  engineering  judgement. 

Specifically,  the  friction  loss  factor  f(L/D)  was  determined 
assuming  f  =  0.003  and  using  the  test  values  of  (L/D) 
provided  in  most  cases.  Reference  26  was  used  to  determine 
the  overall  diffuser  loss  factor  X0. 

As  pointed  out  previously,  it  is  not  easy  to  determine  the 
secondary  entrance  flow  loss  factor  XE.  It  appears  that  the 
intake  geometries  of  the  various  ejector  configurations 
tested  are  of  either  round  or  hellmouth  shrpe,  for  which, 
according  to  the  results  of  Reference  27,  the  entrance  loss 
factor  is  about  0.1.  Thus,  in  the  absence  of  better  infor¬ 
mation  applicable  to  each  specific  configuration,  XE  has  been 
taken  as  0.1  for  the  present  correlation. 

The  correlation  of  the  test  data  with  the  theoretically 
predicted  results  utilizing  the  loss  factors  described  above 
is  performed  in  tabular  form  as  shown  below.  The  reason  for 
the  selection  of  this  method  of  presentation  is  the  fact  that 
the  test  data  obtained  for  a  variety  of  jet  ejector  configura¬ 
tions  are  unsuitable  for  graphical  comparison. 

The  correlations  are  performed  for  the  following  jet  ejector 
configurations : 

A.  SINGLE  CENTRAL  NOZZLE  CONFIGURATIONS 


The  fo1 lowing  references  contain  the  test  data  obtained  for  a 
single  central  nozzle  configuration: 

(i)  Reference  2  presents  the  test  data  for  the  case  of 
an  ejector  without  a  diffuser.  The  correlation  of 
these  data  with  theoretically  predicted  values  is 
given  in  Table  I. 
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TABLE  I 


CORRELATION  OF  THEORY  WITH  TEST  DATA  OF  REFERENCE  2, 
_ SINGLE  NOZZLE  EJECTOR  WITH  NO  DIFFUSER  _ 


aE 

4> 

Tneorv 

4> 

Test 

6.0 

1.30 

1.25 

11.0 

1.36 

1.29 

24. G 

1.42 

1.34 

29.0 

1.44 

1.38 

59.0 

1.51 

1.43 

The  calculated  values  are  based  on  the  following 
loss  factors: 


XE  =  0.1 
f (L/D)  =  0.02 

Xg  ~  0  (No  diffuser) 

The  above  table  indicates  a  satisfactory 
correlation  between  the  theoretical  and  experimental 
results. 

(ii)  The  test  data  obtained  from  Reference  15  are  for  a 
central  nozzle  configuration  without  a  diffuser. 

The  correlation  of  the  test  results  with  the 
theoretically  predicted  values  is  given  in 
Table  II. 
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TABLE  II 


CORRELATION  OF  THEORY  WITH  TEST  DATA  OF  REFERENCE  15, 
S INGLE  NOZZLE  EJECTOR  WITH  NO  DIFFUSER 


aE 

Theory 

4> 

Test 

3.0 

1.22 

1.03 

10.7 

1.36 

1.20 

15.0 

1.37 

1.26 

28.7 

1.44 

1.40 

46.0 

1.48 

1.48 

104.5 

1.55 

1.50 

The  theoretically  predicted  values  are  based  on  the 
following  flow  loss  factors: 

XE  =  0.1 
f (L/D)  =  0.02 

X0  =  0  (No  diffuser) 

The  comparatively  poorer  correlation  at  small  area 
ratios  is  believed  to  be  due  to  the  higher  secondary 
entrance  losses.  Also,  the  test  data  indicate  that 
with  the  increase  of  the  area  ratio  Oe  from  46.0  to 
104.5  (more  than  double),  the  increase  in  augmenta¬ 
tion  ratio  is  only  about  1  percent.  On  the  other 
hand,  from  the  theoretically  predicted  results,  the 
corresponding  increase  in  thrust  augmentation  is 
about  5  peixent.  This  can  be  attributed  to  the 
effect  of  the  nonuniform  velocity  profile  at  the 
secondary  entrance  (not  accounted  for  in  the  theory) 
This  effect  becomes  more  important  at  high  secondary 
to-primary  area  ratios. 


(iii)  Reference  16  quotes  the  following  test  data  for  a 
constant  area  mixing  ejector  without  diffuser: 


a  E 

Test 

129.5 

1.95 

339.0 

2.11 

329.0 

2.21 

369.0 

2.21 

The  thrust  augmentation  ratios  shown  above  are 
either  close  to  or  higher  than  the  maximum  possible 
values  which  can  be  obtained  from  the  idealized, 
one -dimensional  incompressible  analysis,  viz.,  2.0. 
Furthermore,  on  account  of  the  large  area  ratios  of 
aE  >  130,  the  secondary  velocity  at  the  encrance 

plane  would  no  longer  maintain  its  uniformity; 
consequently,  it  would  be  doubtful  whether  the 
maximum  augmentation  ratio  of  2.0  could  be  reached 
at  all  even  if  all  the  losses  were  neglected.  The 
original  work  (Reference  28)  does  not  furnish  any 
usable  information  to  justify  the  unusually  high 
performance  claimed  for  the  simple  central  nozzle 
type  ejector  as  illustrated  therein. 

(iv)  Reference  29  presents  the  test  results  for  the 

ejector  with  and  without  a  diffuser.  The  correla¬ 
tion  of  these  data  with  the  theoretically  predicted 
values  is  given  in  Table  III. 
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TABLE  III 


CORRELATION  OF  THEORY  WITH  TEST  DATA  OF  REFERENCE  29 
SINGLE  NOZZLE  EJECTOR  WITH  AND  WITHOUT  DIFFUSER 

aE 

a  D 

4> 

Theory 

<t> 

Test 

26.36 

1.00 

1.45 

1.37 

7.54 

1.00 

1.34 

1.28 

7.54 

2.49 

1.41 

1.23 

19.55 

1.00 

1.43 

1.35 

5.42 

1.00 

1.29 

1.27 

5.42 

2.49 

1.39 

1.29 

The  test  data  of  Table  III  indicate  that  for 
constant  area  ratios  (  aE=  7.54  or  aE  =  5.42),  the 
addition  of  a  diffuser  with  o0=  2.49  results  in  a 
reduction  of  thrust  augmentation  ratio  <ji>  .  This 
reflects  a  poor  diffuser  efficiency  (high  losses) 
of  the  jet  ejector  configuration  tested. 

(v)  The  tests  of  Reference  30  were  performed  on  a  two- 
dimensional  model  with  no  diffuser.  A  comparison 
of  the  test  data  with  the  theoretically  predicted 
results  is  given  in  Table  IV. 
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TABLE  IV 


CORRELATION  OF  THEORY  WITH  TEST  DATA  OF  REFERENCE  30, 
_ SINGLE  NOZZLE  EJECTOR  WITH  NO  DIFFUSER 


aE 

4> 

Theory 

4> 

Test 

3.24 

1.16 

1.08 

7.13 

1.24 

1.16 

11.5 

1.27 

1.23 

19.0 

1.31 

1.30 

23.0 

1.32 

1.30 

27.0 

1.33 

1.30 

31.0 

1.34 

1.28 

The  predicted  values  are  obtained  by  taking  XE  =  0.1 
and  f(L/D)  =  0.05.  A  higher  friction  loss  factor  is 
used  in  this  case,  as  compared  with  the  above 
configurations,  in  view  of  the  fact  that  due  to  the 
larger  wall  area  in  the  two-dimensional  model,  the 
friction  losses  are  undoubtedly  higher,  especially 
at  low  secondary- to-primary  area  ratios. 

The  results  presented  in  Table  IV  show  a  good 
correlation  between  the  theoretical  and  experi¬ 
mental  values. 

(vi)  Reference  31  presents  the  test  data  on  the  effect 
of  total  primary  stagnation  pressure  ratio  on 
thrust  « ugmentation.  The  test  conditions  are  as 
follows : 
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T0p  =  _>00°K  =  510°R  (50°F) 

Ta  =  283°K  =  478°R  (18°F) 
ag  =  50.0 
L/D  =  5.0 

The  test  data  of  Reference  31  indicate  that  the 
thrust  augmentation  ratio  decreases  gradually 
from  1.24  to  1.22  as  Pop/Pa  increases  from  1.4  to 
1.8,  From  the  analysis,  the  thrust  augmentation 
ratio  is  obtained  as  1.25  by  taking  f(L/D)  =  0.025, 
as  suggested  by  the  paper. 

(vii)  Reference  32  shows  that  at  an  area  ratio  =  2.33, 
the  central  nozzle  ejector  without  diffuser 
(L/D  =  6)  has  its  thrust  augmentation  ratio  <£ 
decreasing  from  1.195  to  1.165  as  the  primary 
pressure  ratio  Pop/Pa  increases  from  1.4  to  1.8. 

The  corresponding  theoretical  value  based  on  loss 
factors  of  f(L/D)  =  0.02  and  Xe  =  0.1  is  1,20. 

B.  MULTIPLE  NOZZLE  CONFIGURATIONS 

The  ♦"est  data  for  the  multiple  nozzle  conf igiirations  are 
obtained  from  the  following  references: 

(i)  The  test  data  reported  in  Reference  13  pertain  to 
an  ejector  configuration  with  multiple  nozzles,  a 
rectangular  mixing  chamber,  and  a  diffuser  with 
lemniscate  contour.  Table  V  shows  the  correlation 
of  the  test  data  with  the  theory. 
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TABLE  V 


CORRELATION  OF  THEORY  WITH  TEST  DATA  OF  REFERENCE  13, 
MULTIPLE  NOZZLE  EJECTOR  WITH  RECTANGULAR 
MIXING  CHAMBER  AND  A  DIFFUSER 


4>  4> 


a  E 

aD 

Theorv 

Test 

35.5 

3.11 

1.38 

1.55 

55.5 

2.73 

1.54 

1.66 

74.4 

2.36 

1.70 

1.79 

135.5 

1.36 

1.79 

1.78 

The  above  table  shows  a  satisfactory  correlation 
between  the  test  results  and  the  corresponding 
theoretical  values.  This  correlation  is  based 
on  the  following  flow  loss  factors: 

X  £  =0.1 

f ( L/D)  =  0.02 

XD  -  0.2 

(ii)  The  test  data  of  Reference  18  and  the  predicted 
values  are  given  in  Table  VI. 


TABLE  VI 


CORRELATION  OF  THEORY  WITH  TEST  DATA  OF  REFERENCE  18, 
_ MULTIPLE  NOZZLES  ARRANGED  IN  A  CIRCLE _ 


a  E 

aD 

Theory 

Test 

14.37 

1.8 

1.71 

1.72 

31.4 

2.2 

1.75 

1.71 

51.0 

2.47 

1.73 

1.90 

The  test  model  used  has  multiple  nozzles  arranged 
in  a  circle.  The  mixing  chamber  is  slightly 
convergent  (the  exit  area  is  8  percent  less  than 
the  entrance  area),  which  does  not  assu  e  constant 
pressure  at  the  walls.  The  predicted  values  shown 
in  Table  VI  are  obtained  using  the  following  loss 
factors : 


Xe  -  0.1 
f (L/D)  =  0.01 

XD  =  0.2 

It  is  seen  that  for  the  first  two  configurations, 
correlation  is  satisfactory,  whereas  for  the  last 
one  the  test  result  is  10  percent  higher  than  that 
theoretically  predicted. 

(ill)  Table  VII  shows  the  effect  of  diffuser  area  ratio 

do  on  thrust  augmentation  for  constant  secondary- 
to-primary  area  ratio  a£  -  12.0.  The  test  data 
were  extracted  from  Reference  21  and  are  applicable 
to  a  four-row  nozzle  configuration. 
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TABLE  VII 


CORRELATION  OF  THEORY  WITH  TEST  DATA  OF  REFERENCE  21, 
FOUR-ROW  NOZZLE  CONFIGURATION  WITH  VARIABLE  DIFFUSER 


aE 

aD 

Theory 

4> 

Test 

12.0 

1.00 

1.37 

1.35 

12.0 

1.19 

1.48 

1.45 

12.0 

1.38 

1.56 

1.55 

12.0 

1.58 

1.59 

1.61 

12.0 

1.77 

1.59 

1.64 

12.0 

2.16 

1.53 

1.57 

Table  VII  shows  a  good  correlation  between  experi¬ 
mental  and  theoretical  results,  indicating  that  the 
best  diffuser  area  ratio  oto  for  optimum  thrust 
augmentation  ratio  is  between  -•  1.5  and  2.0. 
Similar  conclusion  is  reported  in  Section  V,  where 
the  effects  of  a  diffuser  are  discussed. 

(iv)  The  configuration  referred  to  in  Reference  33 
consists  of  a  rectangular  parallel-wall  mixing 
chamber  with  a  diffuser  and  multiple  primary 
nozzles.  The  comparison  of  the  predicted  thrust 
augmentation  ratios  with  the  corresponding 
measured  values  is  presented  in  Table  VIII. 
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TABLE  VIII 


CORRELATION  OF  THEORY  WITH  TEST  DATA  OF  REFERENCE  33, 
SINGLE  AND  THREE-ROW  NOZZLES, 

_ RECTANGULAR  MIXING  CHAMBER  WITH  DIFFUSER 


4> 

Test 

a  E 

aD 

4> 

Theory 

Single-Row 

Nozzle 

Three-Row 

Nozzle 

6.6 

1.23 

1.44 

1.27 

1.36 

12.0 

1.38 

1.56 

1.37 

1.50 

29.3 

1.61 

1.70 

1.45 

1.83 

51.0 

1.73 

1.77 

1.56 

1.89 

The  third  column  in  Table  VIII  represents  predicted 
values  based  on  the  following  values  of  flow  loss 
factors : 


XE  —  0.1 
f (L/D)  =  0.02 
A0  -  0.2 

It  is  seen  that  for  the  first  two  area  ratios 
(aE  =  6.6  and  12,0),  the  predicted  values  are 
higher  than  the  test  data  for  both  configurations, 
viz.,  single  row  of  nozzles  and  three  rows  of 
nozzles.  This  is  believed  to  be  chiefly  due  to  an 
underestimation  of  the  total  head  loss  factor  at 
the  secondary  entrance  X*. 

For  the  other  two  area  ratios  (aE  --  29.3  and  51.0), 
the  predicted  values  are  higher  than  the  test  data 
for  the  single-row  nozzle  configuration  but  lower 
than  those  for  the  three-row  nozzle  configuration. 
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C.  ANNULAR  NOZZLE  CONFIGURATIONS 

Finally,  for  the  annular  nozzle  configuration,  the  following 
test  data  are  utilized: 

(i)  The  tests  as  reported  in  Reference  22  were 

performed  with  single  annular  nozzle  configurations. 

For  the  model  with  a  straight  mixing  chamber 

(no  diffuser),  the  results  are  shown  in  Table  IX. 

TABLE  IX 

CORRELATION  OF  THEORY  WITH  TEST  DATA  OF  REFERENCE  22, 


ANNULAR  NOZZLE  EJECTOR  WITH  NO  DIFFUSER 


aE 

a  D 

.  4> 

Theory 

Test 

13 

1.0 

1.41 

1.29 

17 

1.0 

1.44 

1.38 

22 

1.0 

1.46 

1.42 

The  predicted  values  are  based  on  the  following 
flow  losses: 


X  F  =0.1 

f (L/D)  -  0.01 

\0  0  (No  diffuser) 

The  corresponding  results  for  a  divergent  mixing 
chamber-diffuser  (exit  area  approximate ‘.iy  double 
the  entrance  area)  are  shown  in  Table  X. 
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TABLE  X 


CORRELATION  OF  THEORY  WITH  TEST  DAT 4  OF  REFERENCE  22, 
ANNULAR  NOZZLE  EJECTOR  WITH  DIVERGENT  MIXING  CHAMBER-DIFFUSER 


aE 

5 

Theory 

Test 

8 

1,78 

1.48 

9 

1  v  79 

1.53 

11 

1.81 

1.56 

The  theoretically  predicted  values 

for  the  thrust 

augmentation  ratio  $  have  been  obtained  for  a 
diffuser  area  ratio  of  aD  --  2.0  with  the  following 
loss  factors: 

XE  =  0.1 

f(I./D)  =  0  (No  mixing  chamber) 

X0  =  0.2 

The  large  discrepancy  between  the  actual 
performance  and  the  predicted  results  is  attributed 
to  the  incomplete  mixing  that  can  be  expected  at 
the  exit  of  a  diverging  mixing  chamber-diffuser  of 
relatively  small  length. 

Presented  in  Reference  22  are  also  some  two- 
dimensional  and  three-dimensional  test  results  for 
the  annular  jet  ejectors. 

The  correlation  of  theory  with  the  two-dimensional 
results  is  shown  in  Table  XI. 
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TABLE  XI 


CORRELATION  OF  THEORY  WITH  TWO-DIMENSIONAL  DATA 
OF  REFERENCE  22,  ANNULAR  NOZZLE  EJECTOR 
WITH  DIVERGENT  MIXING  CHAMBER-DIFFUSER 


aE 

aD 

4> 

Theory 

4> 

Test 

9.0 

2.0 

1.80 

1.30 

11.5 

2.0 

1.82 

1.50 

14.0 

2.0 

1.84 

1.43 

The  predicted  values  are  obtained  from  the 
practical  analysis  with  the  following  flow  I^.ss 
factors : 


XE  =  0.1 

f(L/D)  =  0  (No  mixing  chamber) 

X<3  =  0.2 

As  can  be  seen  from  Table  XI,  the  test  results  are 
appreciably  lower  than  the  corresponding  predicted 
values.  This  difference  in  the  results  can  be 
attributed  to  the  following  factors: 

(a)  The  tests  were  pe* formed  for  the  ejector 
configuration  having  a  small  length-to- 
diameter  ratio  (L/D  =3.0).  At  this  ratio 
of  L/D,  it  can  be  expected  that  the 
mixing  at  the  exit  of  the  diffuser  will 
not  be  complete,  contrary  to  the  complete 
mixing  assumption  in  the  theory. 

(b)  Since  the  primary  jet  is  very  close  to  the 
diffuser  walls,  it  is  expected  that  the 
friction  loss  in  the  diffuser  will  be 
increased  as  compared  to  conventional 


central  nozzle  configurations.  This 
increase  in  the  skin  friction  in  the 
diffuser  could  be  accounted  for  by 
increasing  the  diffuser  loss  factor  \D  . 

The  correlation  of  three-dimensional  test  data  of 
Reference  7.2  with  the  theoretically  predicted 
results  is  shown  in  Table  XII. 

TABLE  XII 


CORRELATION  OF  THEORY  WITH  THREE-DIMENSIONAL  DATA 
OF  REFERENCE  22 ,  ANNULAR  NOZZLE  EJECTOR 
WITH  DIVERGENT  MIXING  CHAMBER- DIFFUSER 


a  E 

a  D 

4> 

Theory 

4> 

Test 

20,8 

1.45 

1.66 

1.34 

20.8 

1.63 

1.69 

1.40 

20.8 

1.85 

1.68 

1.47 

20.8 

1.94 

1.66 

1.49 

The  predicted  results  are  based  on  the  same  loss 
factors  as  the  two-dimensional  results.  Also  in 
this  correlation,  the  three-dimensional  test  dat^ 
are  lower  than  the  corresponding  predicted  values 
for  the  same  reasons  as  explained  above  in  the 
two-dimensional  correlation. 

(ii)  Table  XIII  shows  the  correlation  of  the  test  data 

as  presented  in  Reference  34  with  the  theoretically 
predicted  values. 
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TABLE  XIII 


CORRELATION  OF  THEORY  WITH  TEST  DATA  OF  REFERENCE  34, 


THREE-RING  ANNULAR  NOZZLE  EJECTOR  WITH  DIFFUSER 


aE 

aD 

u  ^ 

Theory 

4> 

Test 

44.0 

1.58 

1.93 

44.0 

1.90 

2.02 

2.16 

44.0 

2.30 

2.32 

The  test  data  are  for  the  configuration  with  a 
three-ring  annular  nozzle.  The  theoretically 
predicted  value  of  <£  =  2.02  applies  to  44.0 
and  otp  -  2.0  and  is  based  on  the  following  loss 
factors : 


X^  —  0.1 

f(L/0)  =  0  (No  mixing  euamber) 


Xq  —  0.2 

It  can  be  noted  that  the  theoretically  predicted 
value  is  lower  than  the  test  data  claimed. 
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VII.  RAPID  METHOD  FOR  EJECTOR  PERFORMANCE  PREDICTION 

Presented  in  this  section  is  a  compilation  of  charts  for 
rapid  prediction  of  jet  ejector  performance.  The  numerical 
results  used  in  these  charts  were  obtained  by  solving  the 
theoretical  flow  equations  presented  in  Section  IV. 

A.  IDEALIZED  ANALYSIS 

The  idealized  flow  equations  were  computed  manually  for  both 
constant  area  and  constant  pressure  mixing.  The  corre¬ 
sponding  numerical  results  are  herein  presented  in  Figures 
16  and  17,  which  show  the  variation  of  the  idealized  thrust 
augmentation  ratio  ^  as  a  function  of  secondary- to-primary 
area  ratio  aE  for  a  series  of  constant  values  of  diffuser 
area  ratio  cto* 

B.  PRACTICAL  ANALYSIS 

The  practical  analysis  was  solved  with  the  aid  of  an  IBM 
360  digital  computer  utilizing  FORTRAN  IV  machine  language. 

In  this  case,  two  computer  programs  were  developed,  one  for 
the  incompressible  analysis  including  the  effects  of  major 
flow  losses,  diffuser,  and  forward  speed,  and  the  other  for 
the  compressible  analysis  including  only  the  effects  of  flow 
compressibility.  A  detailed  description  of  the  programs 
including  flow  diagrams  and  typical  computer  outputs  is 
presented  in  Appendix  I.  The  final  computer  results  obtained 
for  static  conditions  are  herein  presented  as  nomographs, 
Figures  18  through  22. 

These  charts  represent  an  effective  analytical  tool  in 
predicting  jet  ejector  performance  including  flow  losses  and 
the  effects  of  flow  compressibility  and  are  suitable  for  use 
in  the  preliminary  design  of  jet  ejectors.  One  of  the 
advantages  of  the  selected  method  of  presentation  is  the 
fact  that  a  wide  range  of  practical  jet  ejector  operating 
conditions  as  well  as  a  variety  of  flow  losses  are  condensed 
in  a  total  of  five  nomographs.  Four  of  these  charts, 

Figures  18  through  21,  contain  the  computer  results  for  the 
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Idealised  Analysis  for  Constant  Pressure 
Mixing. 
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10*0 


Figure  19.  Nomograph  for  Thrust  Augmentation  Ratio  -  Incompressible 
Flow  Analysis  Including  Flow  Losses  (a0  = !  5 , /a  =  0) 


Figure  20.  Nomograph  for  Thrust  Augmentation  Ratio  - 
Flow  Analysis  Including  Flow  Losses  (a0=2.0; 


Flow  Analysis  Including  Flow  Losses  (a0s  2.5 ;  fi  s0) 


incompressible  analysis  (including  flow  losses)  for  the 
diffuser  area  ratios  aj>  =  1.0,  1.5,  2,0,  and  2,5, 
respectively.  The  fifth  nomograph,  Figure  22,  presents  the 
computer  results  for  the  compressible  analysis  for  a  wide 
range  of  the  nondimensionalized  input  parameters  Pa/Pop 
and  Ta/Top 

1.  Evaluation  of  Flow  Losses 

In  utilizing  the  nomographs  for  the  incompressible  analysis, 
the  values  of  the  various  loss  factors  must  be  predetermined. 
No  precise  information  is  available  on  these  loss  factors 
which  are  dependent  on  the  design  and  construction  details 
of  the  system.  However,  the  following  will  serve  as  a 
general  guide: 

a.  Friction  Loss  Factor  f(L/D) 

The  friction  loss  factor  of  a  given  jet  ejector 
configuration  is  a  function  of  the  friction 
factor  f  along  the  ejector  walls  and  the  mixing 
chamber  length-to-diameter  ratio  L/D  required  for 
complete  flow  mixing. 

The  friction  factor  f  =  0.003,  which  is  commonly 
used  for  commercially  smooth  pipes  is  recommended 
to  be  used  for  fairly  smooth  mixing  chamber  walls 
of  jet  ejectors. 

No  precise  information  is  available  for  estimating 
the  mixing  chamber  length-to-diameter  ratios 
required  for  complete  mixing.  However,  in  order 
to  provide  the  designer  some  basis  for  selection 
of  the  required  L/D  ratios,  a  semiempirical 
approach  is  herein  utilized  for  estimating  this 
parameter.  This  approach  is  based  on  the 
assumption  that  for  a  single  nozzle  jet  ejector 
the  flow  losses  due  to  partial  mixing  can  be 
neglected  provided  that  the  mixing  chamber  lei  gth- 
to-diameter  ratio  is  not  less  than  6.0.  This 
implies  that  the  flow  mixing  is  considered  to  be 
complete  for  L/D  >  6.0. 
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The  results  thus  obtained  are  presented  in 
Figure  23.  This  figure  can  be  utilized  to 
estimate  the  total  mixing  chamber  length  required 
for  complete  flow  mixing  for  single,  multiple 
(four  evenly  spaced  nozzles),  and  annular  jet 
ejector  configurations. 

Thus,  utilizing  the  friction  factor  of  f  =  0.003 
and  the  L/D  ratios  from  Figure  23,  the  required 
friction  loss  factor  f(L/D)  for  any  given  jet 
ejector  configuration  can  be  determined. 

b.  Diffuser  Loss  Factor 


The  diffuser  loss  factor  Xq  is  a  function  of  a 
total  head  loss  within  the  diffuser,  the  dynamic 
pressure  at  the  diffuser  entrance  and  its  exit-to- 
entrance  area  ratio.  A  mathematical  definition  of 
this  factor  is  presented  in  the  list  of  symbols. 
Some  usable  data  for  determining  the  loss  factor 
of  various  diffuser  configurations  is  presented  in 
Reference  26. 

c  •  Secondary  Entrance  Loss  Factor 

Very  limited  information  is  available  for 
determining  the  loss  factor  at  the  secondary 
entrance.  This  factor  is  mainly  a  function  of  the 
size  and  shape  of  the  inlet,  but  it  also  depends  on 
the  blockage  effect  of  the  components  of  the 
ejector,  such  as  manifolds,  instrumentation,  etc., 
located  in  the  passage  of  the  secondary  flow 
entrance.  It  is  not  possible  to  detei*mine  this 
factor  accurately,  since  it  varies  from  case  to 
case.  In  Reference  27,  representative  values  of 
"internal  inlet"  loss  factor,  i.e.,  for  static 
conditions  with  no  obstruction  of  the  passage,  are 
presented  as  follows: 
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Figure  23.  Variation  of  Minimum  Mixing  Chamber 
Lengths  Required  for  Complete  Mixing  for  Various 
Ejector  Configurations  (Idealized  Analysis). 
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Type  of  Entrance 


X 


Flared 


Lemniscate  contour 

0.02 

Circular  contour 

0.03 

Straight  contour 

0.05 

Rounded  Edge 

0.1 

Sharp  Edge 

0.3 

It  can  be  noted  from  the  bove  data  that 
lemniscate  lip  contour  i  preferable  in  ejector 
designs.  This  contour  f  ^nds  to  minimize  adverse 
pressure  gradient  ov<  r  ne  lip  and  thus  results  in 
lower  entrance  loss  fa.  .or.  Furthermore s  in 
evaluating  the  entrance  1  ss  factor,  a  consider¬ 
ation  must  be  given  to  the  ratio  of  the  entrance 
diameter  to  the  mixing  tube  diameter.  If  this 
ratio  is  less  than  1.5,  an  additional  loss  in 
thrust  augmentation  ratio  may  occur.  This  loss  may 
be  accounted  for  by  appropriately  increasing  the 
entrance  loss  factor 

2.  Correction  Factor  for  the  Nonuniform  Velocity  Profile  _ 
at  the  Secondary  Entrance 

Section  IV  contains  an  empirical  analysis  for  determining  the 
effect  of  a  nonuniform  velocity  profile  at  the  secondary 
entrance  on  jet  ejector  thrust  augmentation.  This  analysis 
requires  a  knowledge  of  the  flow  nonuniformity  parameter  k 
which  can  only  be  reliably  predicted  from  appropriate  experi¬ 
mental  data. 

The  limited  experimental  data  such  as  present  ?d  in  Table  II 
of  Section  VI  indicates  that  the  approximate  value  of  this 
parameter  is  about  k  -  70.  Thus,  using  Figure  3, the  empirical 
correction  factor  \  for  thrust  augmentation  ratio  due  to 
nonuniform  velocity  profile  at  the  secondary  entrance  can  be 
determined  as  follows: 


(150) 


^K-ro 

X'TT 

where  =  70  is  the  thrust  augmentation  ratio  with  non- 
uniform  secondary  velocity  profile  for  k  =  70,  and  <£x  is  the 
ideal  thrust  augmentation  ratio  for  the  same  aE  . 

3.  Compressibility  Correction  Factor  Cc 

The  compressibility  correction  factor  Cc  can  be  obtained  as 
the  ratio  of  the  compressible  to  the  ideal  value  of  thrust 
augmentation  ratio  <f>  for  the  same  values  of  aE  and  aD  . 
Thus : 


C 


c  = 


(151) 


The  compressible  value  of  thrust  augmentation  ratio  <^>ccan  be 
obtained  from  the  nomograph,  Figure  22,  for  a  given  set  of 
input  conditions  of  Pa /Pop  anc*  •  The  corresponding 

ideal  value  can  be  obtained  from  the  nomograph,  Figure 

18,  using  \E  =  f(L/D)  =  0. 

4  Net  Value  of  Thrust  Augmentation  Ratio 

The  net  value  of  thrust  augmentation  ratio  of  a  given 
ejector  configuration  can  be  expressed  as  follows: 


4>  =XCC^L 


(152) 


where  <*>,  represents  the  incompressible  value  of  thrust 
augmentation  ratio  including  the  effects  of  major  flow 
losses  This  value  can  be  obtained  from  appropriate  nomo¬ 
graphs’,  Figures  18  to  21.  Thus  equation  (152)  can  be  vsed  to 
predict  the  total  thrust  augmentation  ratio  of  a  given 


ejector  configuration  including  the  effects  of  nonuniform 
entrance  velocity,  flow  compressibility,  and  major  flow 
losses  discussed  above. 

5 .  Use  of  the  Nomographs 

The  procedure  for  determining  ejector  thrust  augmentation 
performance  requires  the  use  of  the  nomographs  presented  in 
Figures  18  to  22.  In  order  to  utilize  these  charts, it  is 
first  necessary  to  determine  ejector  geometry,  major  flow 
losses, and  operating  conditions  of  a  given  jet  ejector 
configuration.  The  use  of  the  nomograph  is  explained 
graphically  in  the  key  of  each  chart. 

The  procedure  for  the  use  of  nomographs  for  a  general  case 
cf  an  ejector  with  a  diffuser  is  as  follows: 

(a)  Select  an  appropriate  nomograph  for  a  given 
diffuser  exit-to-entrance  area  ratio  X0. 

(b)  On  the  upper  left-hand  plot  of  the  nomograph, 
draw  a  vertical  line  at  a  given  aE  to  Intersect 
with  the  curve  corresponding  to  the  computed 
value  of  entrance  loss  factor  XE.  Interpolate 
between  XE  curves  if  required. 

(c)  Fi'om  the  point  of  intersection  of  step  (b), 
project  a  horizontal  line  to  intersect  with  the 
curve  corresponding  to  the  computed  value  of  the 
friction  loss  factor  f(L/D).  Interpolate  between 
f(L./D)  curves  if  required. 

(d)  From  the  point  of  intersection  of  step  (c),  draw  a 
vertical  line  to  intersect  the  curve  corresponding 
to  the  computed  value  of  the  diffuser  loss  factor 
X0  .  Interpolate  between  \0  curves  if  required. 

(e)  From  the  point  of  intersection  of  step  (d),  draw 
a  horizontal  line  to  intersect  with  the  scale  of 
the  thrust  augmentation  ratio.  Read  off  this 
value  of  thrust  augmentation  ratio  <£l,  which  in 
this  case  would  correspond  to  the  incompressible 
value  including  the  flow  losses  and  the  effect  of 
a  diffuser. 
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Similar  procedure  is  applied  in  the  use  of  the  nomograph 
for  an  ejector  without  diffuser  (Figure  18),  except  that  the 
vertical  line  of  step  (d)  is  drawn  to  intersect  the 
horizontal  scale  of  thrust  augmentation  ratio.  The  point  of 
intersection  of  the  vertical  line  with  the  scale  yields  the 
required  value  of  <£L  with  no  diffuser. 

The  procedure  for  the  use  of  the  nomograph  (Figure  22)  to 
determine  the  compressible  value  of  thrust  augmentation 
ratio  is  similar  to  that  described  above.  However,  in  this 
case  the  nomograph  is  entered  using  the  precomputed  values 
of  ambient-to-stagnation  pressure  and  temperature  ratios, 
Pa/P0p  and  Ta/T0p  respectively. 

In  the  v  ases  where  the  ejector  performance  is  required  for 
some  intermediate  values  of  diffuser  area  ratios  <*D,the 
usual  interpolation  procedures  between  the  nomographs  can  be 
utilized. 

6 .  Limitations of  the  Nomographs 

Although  the  nomographs  (Figures  18  to  22)  represent  a  rapid 
and  a  practical  analytical  tool  in  evaluating  performance 
of  a  given  jet  ejector  configuration,  the  usefulness  of  the 
charts  is  limited  by  the  assumptions  inherent  in  the 
analysis. 

lire  major  problem  exists  in  the  user's  ability  to  accurately 
predetermine  the  required  flow  losses  of  a  given  ejector 
configuration. 

The  assumptions  used  to  determine  the  friction  factor  along 
the  mixing  chamber  walls  (f-  0.003),  the  total  mixing  chamber 
length  for  complete  mixing, and  the  flow  nonuniformity 
parameter  (  k-  70)  for  the  velocity  profile  at  the  secondary 
entrance  require  further  experimental  verification 

The  nomographs  have  been  carefully  prepared  and  their 
accuracy  is  expected  to  be  within  percent  of  the 
computer  results. 


C .  PROCEDURE  FOR  DETERMINING  PERFORMANCE  OF  /  GIVEN  JET 
EJECTOR  CONFIGURATION 


The  following  procedure  can  be  utilized  to  determine  the 
performance  of  a  given  jet  ejector  configuration: 

(i)  Determine  the  following  geometric  parameters 
and  operation  conditions: 

(a)  Type  of  ejector  configuration 

(b)  Number  of  primary  nozzles,  N 

(c)  Secondary- to-primary  area  ratio, 

(d)  Diffuser  exit- to-entrance  area  ratio,  aD 

(e)  Diffuser  length  and  expansion  angle 

(f)  Mixing  chamber  shape 

(g)  Ejector  intake  geometry 

(h)  Ejector  operating  conditions,  Pg/Pop 
and  Ta/T0p 

(ii)  Knowing  and  the  type  of  ejector  configuration, 
enter  Figure  23  and  obtain  the  total  mixing 
chamber  length  (L/D)  required  for  cc—Tete  mixing. 
For  the  case  of  a  multiple  nozzle  ejector  configu¬ 
ration  with  N*  4.0,  determine  the  total  mixing 
chamber  length  using  the  following  equation: 


t— )=  -4=(— )  053) 

D  m  vN  D  § 


where  (L/D)s  can  be  obtained  from  Figure  23. 


(iii)  Assuming  friction  factor  f  =  0.003,  compute 
ejector  friction  loss  factor  f(L/D)  using 
(L/D)  value  from  step  (ii). 

(iv)  Using  data  of  Reference  26  (or  other  pertinent 
data)  and  the  diffuser  geometry  from  step  (i), 
determine  diffuser  loss  factor  XD  . 

(v)  Using  data  of  Reference  27  (or  other  pertinent 
data)  and  ejector  intake  geometry  from  step  (i), 
determine  the  ejector  entrance  loss  factor  XE  . 

(vi)  With  the  flow  losses  determined  in  steps  (iii)  to 
(v)  and  known  values  of  aE  and  aD  from  step  (i), 
enter  the  appropriate  nomograph  and  obtain  the 
incompressible  value  of  thrust  augmentation  ratio 
including  flow  losses,  <pL  . 

(vii)  Using  aE  from  step  (i)  an'  assuming  k  =  70,  enter 
Figure  3  and  obtain  (pK.7Q  and  <pz  .  Then 
compute  the  empirical  correction  factor  x  from 
equation  (150).  This  factor  accounts  for  the 
reduction  of  thrust  augmentation  ratio  due  to 
nonuniform  velocity  profile  at  the  secondary 
entrance. 

( viii)  Us'-.-v*  the  ' Uji tio.  Oi.  Pq ^ ^op  a*tu 

and  aE  from  step  (i),  enter  nomograph  Figure  22 
and  determine  the  compressible  value  of  thrust 
augmentation  ratio  cf>Q  .  Also,  enter  Figure  18 
using  X E  -  f(L/D)  -  0  and  obtain  ideal  value  of 
thrust  augmentation  ratio  .  Then  compute  the 
compressibility  correction  factor  Cc  using 
equation  (151). 

(ix)  Using  <£l  from  step  (vi),  x  from  step  (vii),  and 
Cc  from  step  (viii),  compute  the  required  thrust 
augmentation  ratio  from  equation  (152). 

(x)  For  a  known  mixing  chamber  shape, compute  mass 
entrainment  ratio  vr  using  equation  (34)  for 
constant  area  mixing  or  equation  (51)  for  constant 


1 15 


pressure  mixing.  The  mass  entrainment  ratio  for 
compressible  analysis  or  incompressible  analysis 
including  major  flow  losses  can  be  most  conven¬ 
iently  obtained  using  the  computer  program 
described  in  Appendix  I. 

D.  SAMPLE  CALCULATION 

To  more  clearly  indicate  the  analytical  procedures  for 
determining  ejector  performance,  a  sample  calculation  is 
performed  as  follows: 

(i)  Assume  a  single,  central  nozzle  configuration 

of  constant  area  mixing  chamber  having  aE  =  100 
and  a0  =  2.0.  Also  assume  the  following 
operation  conditions  Ta  =  70°F,  Tqp  =  600°F, 

Pa  --  14.7  Ib/in^and  P0p  -  21  lb/in^. 

Compute 


To  .  70+460  .  Jj30 
T0~' 600+460  '1060 


Pq  14.7 
P0p  '  21.0 


0.7 


(ii)  Using  aE  =  100, enter  Figure  23  and  obtain 


(-5-)=  6.43 
u  s 

(iii)  Assuming  friction  factor  f  -  0.003,  calculate 

ejector  friction  loss  factor  using  (L/D)  -  6.43 

from  step  (if)*  Thus, 


f (L/D)  =  0.003  X6.43  =  0.0I93 
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(iv)  Using  data  of  Reference  26,  determine  the 
diffuser  loss  factor 


Xd=0.2 


(v)  Using  data  of  Reference  27,  determine  entrance 
loss  factor 


XE=  0  1 


(vi)  With  the  flow  losses  computed  in  steps  (iii)  to 
(v)  and  using  aE  =  100,  atp  =  2.0, enter 
nomograph  Figure  20  and  determine  the  incompressible 
value  of  thrust  augmentation  ratio  including  flow 
losses 


<pLz  1.98 


vvri)  Assuming  k  ~  70  and  using  aE  =  100,  enter 
Figure  3  and  obtain 


^K=70= 


i.745 


<^>x=  1.776 


Compute  the  empirical  correction  factor  %  for 
thrust  augmentation  ratio  due  to  nonuniform 
velocity  profile  at  the  secondary  entrance;  thus, 
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( viii) 


1.745 

1776 


=  0.983 


Using  the  values  Rj/Pop  0.7  and  To/ Top  =  0*5 
computed  in  step  (i),  enter  nomograph  Figure  22 
and  obtain  the  compressible  value  of  thrust 
augmentation  ratio  corresponding  to  aE  =  100. 
Thus , 


<f>c-  1.737 


Also  enter  nomograph  Figure  18  using 

XE  =  f (L/D)  =  0  and  aE  =  100  and  determine 

ideal  value  of  thrust  augmentation  ratio.  Thus, 


<£,=  1.776 


Compute  compressibility  correction  tactor 


Cc  = 


h. 

4>i 


1.737 

1.776 


0.978 


(ix)  Finally,  using  <£u  =  1.98  (from  step  (vii), 

X  =  0.983  (from  step  (vii)),  and  Cc  =  0.978 
(from  step  viii),  compute  the  required  thrust 
augmentation  ratio  as  follows: 


4> :  4>l  Xcc  =  198  *  0  983  x  0  978  5  190 
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(x)  Using  ae  =  100  and  do  =  2.0,  compute  mass 

entrainment  ratio  w  from  equation  (34);  thus, 


•198  +  100. 


:I0Q04 


10004 


15.49 
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VIII.  CONCLUSIONS  AND  RECOMMENDATIONS 


1.  A  review  of  the  technical  literature  shows  that  the 
existing  analytical  and  experimental  information  on 
the  thrust  augmentation  characteristics  of  jet  ejectors 
cannot  be  used  as  an  effective  design  tool. 

2.  The  analysis  performed  in  this  report  indicates  that  a 
constant  area  mixing  ejector  yields  a  higher  thrust 
augmentation  ratio  than  an  equivalent  constant  pressure 
mixing  configu  ation. 

3.  For  any  practical  value  of  secondary- to-primary  area 
ratio,  the  thrust  augmentation  ratio  reaches  an  optimum 
value  with  the  diffuser  area  ratio  ranging  between  1.5 
and  2.0. 

4.  To  ootain  a  maximum  tnrust  augmentation,  the  mixing 
chamber  length  should  be  compromised  so  as  to  achieve 
the  best  mixing  with  a  minimum  wall  friction. 

5.  Annular  and  multiple  nozzle  ejectors  require 
substantially  shorter  mixing  chamber  lengths  for  complete 
mixing  as  compared  to  an  equivalent  single,  central 
nozzle  configuration. 

6.  Tb.  fiw  1  >s.  es  have  predominant  effect  on  jet  ejector 
pei  .-ox manee ,  whereas  flow  compressibility  is  only  of 
secondary  importance. 

7.  An  increase  of  forward  speed  (parallel  or  perpendicular 
to  the  ejector)  causes  a  decrease  of  the  thrust 

aut>  lentation  ratio. 

8.  The  available  test  data  are  insufficient  to  determine 
reliably  the  validity  of  the  assumptions  utilized  in 
the  analyses.  It  is  therefore  recommended  that  a 
systematic  test  program  be  conducted  to  determine  the 
applicability  of  these  assumptions  as  well  as  to 
provide  more  precise  information  to  evaluate  the 
empirical  correction  factors  presented  in  this  report. 
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APPENDIX  I 


COMPUTER  PROGRAMS  FOR  SOLUTION  OF  FLOW  EQUATIONS 
FOR  THE  PRACTICAL  ANALYSIS 


Due  to  the  complexity  of  the  flow  equations  which  require 
tedious  and  lengthy  computation,  it  was  decided  that  a  more 
efficient  and  accurate  method  of  obtaining  the  required 
results  would  be  by  the  use  of  a  digital  computer.  The 
practical  analysis  presented  in  Section  IV  was  therefore 
programmed  on  an  IBM  360  computer  utilizing  BPS  FORTRAN  IV 
language.  Two  separate  prog. ..ms  were  performed.  The  first 
program  deals  with  the  incompressible  flow  analysis  whereby 
major  flow  losses  and  the  effect  of  a  diffuser  and  forward 
speed  are  included.  The  second  program  deals  exclusively 
with  flow  compressibility. 

1 .  Incompressible  Flow  Analysis 

As  mentioned  previously,  the  effect  of  the  major  flow  losses, 
diffuser,  and  forward  speed  are  investigated  on  the  basis  of 
the  incompressible  flow  analysis.  This  is  accomplished  by 
solving  equation  (57)  for  the  velocity  ratio  V^/V,  .  The 

velocity  ratio  V|S/Vlj3  is  then  computed  from  equation  (58). 
The  thrust  augmentation  ratio  and  mass  entrainment  ratio  w 
are  obtained  from  equations  (60)  and  (62),  respectively.  A 
simplified  flow  diagram  of  this  program  is  presented  in 
Figure  24,  and  a  typical  IBM  computer  output  is  shown  in 
Table  XIV. 

2.  Compressible  Flow  Analysis 

The  compressible  flow  analysis  as  presented  in  Section  IV 
involves  a  simultaneous  solution  o.  the  two  nonlinear 
equations  (85)  and  (86)  for  the  two  unknowns  V2/VQ  and 
Tjp/Tap  in  terms  of  given  input  parameters  Pa/P0p >  'a/T0p  , 
and  Oe*  However,  in  generating  the  required  computer  data, 
equations  (85)  and  (86)  are  solved  in  their  equivalent 
dimensional  form  as  follows: 


125 


Figure  24.  Computer  Flow  Diagram 
for  Incompressible  Analysis. 
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l%^Pd-ftpCpv[(a;+!>^(^^)  1-0 


...(154) 
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The  two  nonlinear  equations,  (154)  and  (155),  are  solved  j.or 
the  two  unknowns  and  V2  in  terms  of  dimensional  input 
values  of  Pq,  Pop  ,  %  ■  and  T0p.  The  solution  of  these  two 
equations,  designated  as  if/,  and  \pz>  respectively,  is 
obtained  using  the  Newton-Raphson  iteration  procedure  and 
Taylor’s  series  as  follows'. 


(Vi+.=  Vi  +th)i 


(156) 


lV2)i  +  ,  =  (V2)i  +(k)i 


(157) 


where  i  =  0,  1,  2  ...  representing  successive  iterations. 
Using  Taylor’s  series,  the  functions  (^()j+|  and 

)  can  be  expressed  as  follows: 


‘t'-WW.  +lh'i  (ad  +w.b* 


(158) 


/d%\  (d'l'z 

+(kfe 


(159) 


Solving  for  (h);  and  00 j  from  equations  (158)  and  (159) 
yields 

df, 


dVj  ; 

(h>r  dji^ 

dp{p 

d^2  d\f/z 


(160) 


tylD  dV2  « 


i^'P 
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di>, 

a^  "*■ 

d*2 


1  it  it 

dp>  p  dv2 

it  it 

<3/>ip  av2 


(161) 


The  above  i  teration  procedure  is  started  assuming  initial 
conditions  (i  r-  0) 


lVo  =  K'^P 

(162) 

<V2>„  =  K2VQp 

(163) 

where  K,  and  K2  are  suitable  programmed  constants  and  pa p 
and  Vap  are  obtained  from  equations  (71)  and  (72), 
respectively. 

A  simplified  flow  diagram  of  this  computer  program  is  shown 
in  Figure  25,  and  a  typical  sample  of  final  output  is 
presented  in  Table  XV. 
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